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ABSTRACT 

Based on high-resolution spectra obtained during gravitational microlensing events we present a detailed elemental abundance analysis 
of 32 dwarf and subgiant stars in the Galactic bulge. Combined with the sample of 26 stars from the previous papers in this series, 
we now have 58 microlensed bulge dwarfs and subgiants that have been homogeneously analysed. The main characteristics of the 
sample and the findings that can be drawn are: (i) The metallicity distribution (MDF) is wide and spans all metallicities between 
[Fe/H] = -1.9 to +0.6; (ii) The dip in the MDF around solar metallicity that was apparent in our previous analysis of a smaller 
sample (26 microlensed stars) is no longer evident; instead it has a complex structure and indications of multiple components are 
starting to emerge. A tentative interpretation is that there could be different stellar populations at interplay, each with a different scale 
height: the thin disk, the thick disk, and a bar population; (iii) The stars with [Fe/H] < -0. 1 are old with ages between 10 and 12 Gyr; 
(iv) The metal-rich stars with [Fe/H] > -0.1 show a wide variety of ages, ranging from 2 to 12 Gyr with a distribution that has a 
dominant peak around 4-5 Gyr and a tail towards higher ages; (v) There are indications in the [a/Fe] - [Fe/H] that the "knee" occurs 
around [Fe/H] = -0.3 to -0.2, which is a slightly higher metallicity as compared to the "knee" for the local thick disk. This suggests 
that the chemical enrichment of the metal-poor bulge has been somewhat faster than what is observed for the local thick disk. The 
results from the microlensed bulge dwarf stars in combination with other findings in the literature, in particular the evidence that 
the bulge has cylindrical rotation, indicate that the Milky Way could be an almost pure disk galaxy. The bulge would then just be a 
conglomerate of the other Galactic stellar populations (thin disk, thick disk, halo, and ...?), residing together in the central parts of the 
Galaxy, influenced by the Galactic bar. 
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1. Introduction 

Galactic bulges are emerging as inherently complex features in 
spiral galaxies an d are a defining component of many spiral 
and d isk galaxies dKormendv & Kennicutjl2004l : lAthanassoulal 
I2005h . Numerous studies have shown them to have several spa- 
tial structures overlaying each other. The Milky Way is no dif- 
ferent and over the last decades studies of the Galactic bulge 
have revealed an increasingly complex, and not always consis- 
tent, picture of the stellar population present in the innermost 
few kpc of our Galaxy. The formation of this central concen- 
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tration of stars is not fully understood and if the observed spa- 
tial structures are uniquely related to dynamical and chemical 
features is a very actively studied field. New surveys such as 
the VISTA Varia bles in the Via Lactea (VVV) public survey 
dSaito et al.ll2010l) are looking into this intriguing Galactic com- 
ponent. Several outstanding questions remain un-resolved in- 
cluding the true shape of the metallicity distribution function 
(MDF) and how the MDF is connected to different spatial and 
kinematical stru ctures. A recent discuss ion of these issues can 
be found in, e.g jBabusiaux etaDd2010h . 

A number of formation scenarios for the G alactic bulge have 
been p roposed and were first summarised by IWvse & Gilmorel 
( 1992b and set out in grea t er detail in later w orks (e.g., 
iKormendv & Kennicuttl 120041: lAthanassoulal 120051) . The first 
major scenario is in situ star formatio n from primordial 
gas in the gravitational collapse (e.g., lEggen et aD Il962t 
iMatteucci & Brocatolll990j). or merging of clumps within th e 
disk at high redshift (e.g.. lNoguchill999blBournaud et al.ll2009l) . 
Historically, this has been thought to be the major component 
of the Galactic bulge, and bulges in general, and is referred to 
as a classical bulge. However, in the Milky Way it now seems 
that the spherical component that should form from an initial 
collapse is after all a rather small component. The mass of the 
classical bulge m ight be less than 8 % of the mass of the disk 
(IShen et al.ll20loh . 

Instead, the scenario that is currently emerging as the most 
compelling is secular evolution in which a bulge and bar forms 
from buckling instabilities in the disk (e.g., Combes et al.lll990t 
Kormendv & Kennicuttl l2004t lAthanassoulal 120051: iNess et all 
20121) . Bars are common in other galaxies and the exis- 
tence of a bar in the Milky Way in addition to its classical, 
spherical bulge has been debated for som e time. Star counts 
have shown that there definitely is a bar ([Stanek et"aT] 1 1 994t 
iMcWilliam & Zoccalill2010l:ISaito et al"1l201 lll2012l) . Studies of 
the stellar kinematics show tha t a large fraction of the stars ex- 
hibit cylindrical rotatio n (e.g., IShen et all 120101 : iHoward et al.l 
l2009HSumi et"aT1l2003h which is ind icative of a boxy bulge, i.e., 
one with a bar ( lAthanassoulal 120051) . Thus both star counts as 
well as the stellar kinematics unambiguously show the presence 
of a bar. Since the bar forms from buckling instabilities in a disk, 
there may be a strong connection between the Galactic stellar 
disk and the bar population. 

A realistic approach could be a mixed scenario encom- 
passing two or more formation scenar ios. Mixed formation has 
been observed in other galaxies (e.g. JPrugniel et al1l2001l) and 
a mixed origin of the Galactic bulge is explored by a number 
of theoretical studies predicting both the shape of the MDF as 
well as the elemental abundance trends that should be expected 
(e.g.. iTsuiimoto & Bekkill2012t iGrieco et al.ll2012l) . The pres- 
ence of the first s t ars in the Galactic bulge as predicted in, e.g., 
IWvse & Gilmorel l[l992) who discuss the possibility of merging 
of dwarf satellite galaxies that sink into the Galactic bulge, have 
been furt her developed wi thin the context of cosmological sim- 
ulations jTumlinsonl20 1 Oh . 

Elemental abundances in stars are able to provide vi- 
tal clues to the formation time scales involved in the stel- 
lar population(s) present in the Galactic bulge as well as 
to c onstrain the shape of the MDF and slope of the IMF 
(e.g jBalrero et al]l2007l ITsuiimoto & Bekkil2012l) . Initial stud- 
ies of red giants revealed the Galactic bulge to be rela- 
tively metal-rich but with a wide metallicity distribution (e.g., 
IMcWilliam & Richlfl99l IZoccali et al1l2008l) . It also appeared 
enhanced in the q-elements, indicative of a short forma tion time 
(IMcWilliam & Richlll994t iFulbright et alll2006l 120071) . The el- 



evated levels of a-elements continued to super-solar metallici- 
ties, thus setting the bulge stellar population apart from what was 
found in the solar neighbourhood, i .e., with solar levels of m ost 
elements at these metallicities (e.g jEdvardsson et al.ll 1 993h . As 
the bulge is far away the stars are faint and hence these stud- 
ies were and, essentially are, confined to the most luminous gi- 
ants. Results based on giant spectra are not trivial to interpret. 
Physical processes within giant stars also erase some original 
abundance signatures. This is the case for C, N, and Li in all gi- 
ants, and for O, Na, Mg, and Al in some cases (e.g.. lKraft et al.l 
|1992|) . Additionally, although recent differential analysis of gi- 
ant stars in the bulge and disk have achieved high precision 
dAl ves-Br ito et al .1201 Oh . the situation at high metallicities is still 
unclear due to heavy blending, as the cool atmos pheres of giants, 
rich i n molecules, are difficult to analyse (e.g., iLebzelter et al.l 
120121) . 

The spectra of dwarfs, even metal-rich ones, are fairly 
straightforward to a nalyse and are the best tracers of galactic 
chemical evolution (lEdvardsson et al.l Il993l) . The dwarf stars 
are also unique in that they give star-by-star age information 
and not just the bulk age-information available through analy- 
sis of colour magnitude diagrams. Normally dwarf stars in the 
bulge are too faint to be observed with high-resolution spectro- 
graphs. However, during microlensing events they can brighten 
by several magnitudes and during the last few years there 
have been several studies that have utilised this phenomenon 
and obtained spectra that allowed detailed elemental abundance 
analysis of bulge dw arf stars. In total t here have been 26 
bulge dwarfs observed (|Cavallo et al.ll2003 : Johnson et al 1I2007L 



2008l | Cohen et al.ll2008l 120091: iBensbv et al.l2009HEpstein et al.l 
2010HBensby et al. | 201 Obi 1201 ll) all homogeneously analysed in 
Bensbv et al.l ( 201 ll) . These studies have revitalised the study of 



chemical evolution and structure of the bulge. 

For instan ce, the microlensed dwarf stars dBensbv et al.l 
1201 Obi I201T1) as well as of red gian t branch (RGB) stars 
dHill et all 1201 ll lUttenfhaler et all 1201 2|) have revealed a dual- 
component bulge MDF - one metal-poor component around 
[Fe/H] as -0.6 and one metal -rich component around [Fe/H] « 
+0.3. There has also been considerable evidence that secular 
evolution of the thick disk may have played a large role in the 
origin and evolution of the bulge. This is based on the apparent 
similarity of the [a/Fe] - [Fe /H] abundance trends in the bulge at 
sub-so l ar metallicities (e.g..lMelendez et al. 2008; Be nsbv et al.l 
l2010bl lAlves-Brito et al1l2010l iGonzalez et al.ll201 ll) and those 
of the nearby thick disk (e.g, Bensb vet al.l teOOTaTlFuhrmannl 
120081 iReddv et aDl2006l) . Also, the metal-poor peak in the bi- 
modal bulge MDF appears to coincide with the MDF for the 
thick disk and considering that the age structure of the metal- 
poor bulge dwarfs is the same as for the thick disk, i.e., mainly 
old, it might well be that the bulge and t he thick disk have had 
similar or even shared chemical histories dMelendez et al.ll2008l 
Alves-Brito eTaT] 120101 IBensbv et alJ 1201 Obi Gonzalez et al. 



201 ll IBensbv et al.ll201ll) . Of course, the comparisons between 



the bulge and the thick disk have mainly used nearby thick disk 
samples, and if the bulge partly originated through secular evolu- 
tion of the thick disk it would be desirable to compare with thick 
disk stellar samples locat ed in situ in the inne r Galactic disk. 
The only such study is the IBensbv et al.l d2010al) study of 44 red 
giants, located between galactocentric radii 4-7 kpc and up to 
3-4 kpc from the Galactic plane. Those results show that the in- 
ner disk region has a similar abundance structure as seen in the 
solar neighbourhood, i.e., a thin disk-thick disk dichotomy, and 
that the similarities between the bulge and the thick disk persist 
when comparing to an inner disk sample. 
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Table 1. Summary of the 32 microlensed bulge dwarf stars new in this study (sorted by observation date) 1 . 
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Notes. ' Given for each microlensing event is: RA and DE coordinates (J2000) read from the fits headers of the spectra (the direction where the telescope pointed 
during observation); galactic coordinates (/ and b); duration of the event in days (Te); time when maximum magnification occured (T max ); maximum magnification 
(A„); time when the event was observed with high-resolution spectrograph (T B i, s ); the exposure time (Exp.), the measured signal-to-noise ratio per pixel at 
~6400 A (S/N); the spectrograph that was used (Spec); the spectral resolution (/?). 



A recent controversy is the overall age distribution in the 
bulge. Photometric observations of different bulge fields point 
to an e xclusively old stellar population with ages greater than 
lOGyr JZoccali etal l l2003t IClarkson et all 120081 iBrown et al.l 
l2010t IClarkson et al.l l201lh . In addition to the very old stars, 
the observations of microlensed bulge dwarf stars have re- 
vealed a significant fraction of stars with ages betw een 3-7 Gyr 
dBensbv et al.ll2010bl 120111) . iNataf & Goulrj (1201 2|) tried to rec- 
oncile the dilemma of the young bulge dwarfs by showing that 
a factor of 2 discrepancy between spectroscopic and photo- 
metric age determinations of the Galactic bulge main-sequence 
turnoff possibly can be explained if the Galactic bulge is helium- 
enhanced relative to that assumed by standard isochrones. An 
intriguing recent result, based on spectra of 575 bulges from the 
Sloan Digital Sky Survey, is that barred galaxies show a stellar 
age distribution w ith a young stellar comp onent not present in 
unbarred galaxies dCoelho & Gadottill201 lb . 

Another intriguing result is the multi-component (> 4) MDF 
proposed by Ness et al. (in prep.) based on low/intermediate- 
resolution spectra of 28 000 RGB stars at different longitudes 
and latitudes in the bulge. Whether these proposed multiple 
bulge components are real and where they come from needs to 
be investigated. 



In this work we present a detailed elemental abundance anal- 
ysis of 32 new microlensing events towards the Galactic bulge. 
T ogether with the prev ious 26 events homogeneously analysed 
in lBensbv et al.l (1201 ll) . the sample now contains 58 microlensed 
dwarf and subgiant stars in the bulge. 

2. Observations and data reduction 

Dwarf and main sequence turn-off stars i n the Galactic bulge 
have V magnitudes around 19-20 (e.g., iFeltzing & Gilmorel 
2000), which is far too faint to observe under normal circum- 
stances even with today's 8-10m class telescopes. To obtain a 
decent high-resolution spectrum (S /N > 50 and R > 40 000) of 
V — 19-20 star would require more than 50 hours of exposure 
time with an instrument such as UVES on VLT. To the rescue, 
gravitational microlensing has surfaced as nature's golden mag- 
nifying glass to study dwarf stars in the bulge. Gravitational mi- 
crolensing occurs whenever a compact object, by chance, passes 
the line of sight between an observer and a distant object. The 
brightness of the distant object can then be magnified by factors 
of several hundred. However, as microlensing events are com- 
pletely random it is impossible to predict when and where they 
will happen. The search for microlensing events in the bulge 
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Fig. 1. Light curves for the first 16 of the 32 new microlensing events. Each plot has a zoom window, showing the time intervals 
when the source stars were observed with high-resolution spectrographs. In each plot the un-lensed magnitude of the source star is 
also given (7$). 



is today mainly done by the MOA and OGLE surveys, which 
monitor the bulge every night. The OGLE survey is carried 
out primarily in the / band using 1.3 m Warsaw Telescope at 
Las Campanas, Chil e and detects ro ughly 1700 events per year 
dUdalskiet al.l Il994t lUdalskil 12003b . The MOA survey is car- 
ried out primarily in a broad R/I filter using the 1.8 m telescope 
at M t. John, New Zealand, and detects r oughly 650 events per 
year jBond e t al. 200ll lSumi et"aIll20Tll) . For a small subset of 



these events the background source stars are both (1) dwarfs (or 
subgiants), and (2) magnified highly enough to allow us to ob- 
tain a high-resolution high-S/N spectrum in a few hours. To find 
these spectroscopy targets, we perform real-time modelling of 
OGLE and MOA data as they are updated, and also acquire ad- 
ditional data on promising targets using the ANDICAM camera 
on the 1.3 m SMARTS telescope at Cerro Tololo Interamerican 
Observatory (CTIO) in Chile. We also use the SMARTS tele- 
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Fig. 2. Light curves for the remaining 16 of the 32 new microlensing events. Each plot has a zoom window, showing the time 
intervals when the source stars were observed with high-resolution spectrographs. In each plot the un-lensed magnitude of the 
source star is also given (Is,). 



scope to measure the (V - I) colours of most events when they 
are near peak brightness. 

Due to the unpredictability of microlensing events it is im- 
possible to have a standard observing program, which is allo- 
cated specific nights in visitor mode or put in queue service 
mode. Hence, starting in the ESO observing period P82 we 
have be en running Target-o f-Opportunity (ToO) programs using 
UVES (iDekker et al.l2000h on the VLT allowing us to trigger ob- 



servations with only a few hours notice. Up until now, including 
the second of five triggers in P90, we have submitted 52 trig- 
gers with the result of 47 successfully acquired spectra of dwarf 
stars in the bulge. Of the five "failures" there were two events for 
which the predicted magnification was much higher than the ac- 
tual one, resulting in too low S /N and spectra useless for abun- 
dance analysis, one event turned out to be a red giant and will 
be published in another paper, one spectrum was obtained when 
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only the UVES blue arm was available resulting in very low S /N 
and short wavelength coverage, and deemed use less, and one 
event turned out to be a nova towards the bulge (see lWagner et alj 
120121) . The ToO observations with VLT have been augmented 
with a few observatio ns obtained with the MIKE spectrograph 
( Bernstein et alj 2003b on Magellan or the HIRES spectrograph 
( Vogt et al.ll994 ) on Keck. The spectra have resolutions between 
R « 40 000 - 90 000 (see Table[T}. 

Our goal has been to observe the targets as close to peak 
brightness as possible, and to observe targets that reach at least 
/ « 15.0 during the observations. With a 1" wide slit this means 
that the spectra should have a minimum signal-to-noise ratio 
of around S /N ~ 50 - 60 and a resolution of R ~ 42 000 
in the UVES 760 nm setting (according to the UVES exposure 
time calculator). However, sometimes the observed microlensing 
light curve does not follow the predicted one, meaning that the 
obtained spectra could be better or worse than predicted. Also, 
observing conditions (seeing and clouds) have severely affected 
the quality of some of the observations. Hence, the signal-to- 
noise ratios vary between S /N ~ 30 - 100 (see Table [TJ and for 
a few events the S/N is even lower than 20. Figure Q] shows the 
light curves for the 3 2 new events (light curves for the first 26 
events are published in lBensbv et al.l2009ll20To5l201 lb . In each 
plot the time interval when the spectroscopic observations were 
carried out has been marked (see insets in each panel). Generally, 
we were able to catch many objects very close to peak bright- 
ness, but in a few cases we were less lucky. Especially interest- 
ing is the observation of MOA-2012-BLG-532S (Fig. 13. This 
event suddenly entered a caustic crossing resulting in a strong 
a sudden increase in brightness. We managed to submit a trig- 
ger and start observations. However, its brightness started to de- 
crease again sooner than expected, resulting in most of the two 
hour exposure being executed when it had dropped 2 magnitudes 
in brightness. Hence, the spectrum for MOA-2012-BLG-532S is 
of lower quality. 

Data reductions of the UVES spectr a were carried out with 
the UVES pipeline dBallester et al.ll2000t) versions 4.9.0 through 
5.0.7, depending on when the event was observed. The MIKE 
spectra were reduced with the Carnegie Observatories python 
pipelineQ and the red uction of the HIRE S spectrum follows the 
procedure outlined in lCohen et alj (120081) . 



3. Analysis 

3. 1 . Stellar parameters - with Fe I NLTE corrections 

NLTE corrections for Fe i are now included in the analysis com- 
pared to the anal ysis of the pre vious s ample of 26 microlensed 
dwarf stars from lBensbv et all fcOlObl |2011|) . As the inclusion 
of NLTE corrections for Fe i lines impacts the final values of 
the stellar parameters, those 26 stars have been re-analysed, en- 
suring a homogeneous treatment of the full sample. The NLTE 
corrections are implemented on a line-by-li ne basis using a n 
IDL script kindly provided by Karin Lind dLind et alj|2012l) . 
Fen lines are on the other hand not sensitive to NLTE effects 
(e.g jThevenin & Idiarlll 1 999b iLind et all2012l) . and need not be 
corrected. Apart from the implementation of Fe i NLTE correc- 
tions, the methodology to find the s tellar parameters is the same 
as be fore and is fully described in iBensbv et al.l d2009l 1201 Obi 
1201 lb . In brief, the analysis is based on sta ndard 1-D plane- 
parallel MARCS model stell a r atmospheres ( Gustafsson et al.l 
1 19751 lEdvardsson et al] 119931: lAsplund et al.lll997h . and abun- 
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Fig. 3. Comparison of stellar parameters before and after imple- 
mentation of NLTE corrections for Fe i lines. The differences are 
NLTE values minus LTE values and mean differences and <x are 
indicated in each panel. 



dances are calculated with the Uppsala EQWIDTH program us- 
ing equivalent widths measured by hand using the IRAF task 
SPLOT. The effective temperatures were determined from exci- 
tation balance of NLTE corrected abundances from Fe i lines, 
surface gravities from ionisation balance between NLTE cor- 
rected Fe i abundances and abundances from Fe n lines, and the 
microturbulence parameters by requiring that the NLTE cor- 
rected Fe i abundances a re independent of line strength. 

As demonstrated bv lLind et ail d2012l) the NLTE effects are 
very small, being smallest for the most metal-rich stars (gener- 
ally below 0.01 dex) and increasing for lower metallicities, and 
for stars with lower surface gravities. Figure [3] shows compar- 
isons of effective temperatures, metallicities, and microturbu- 
lence velocities for all 58 stars with and without the inclusion 
of Fe i NLTE corrections, and as can be seen the effects are truly 
minuscule for the type of stars analysed here (mostly metal-rich 
F and G dwarfs). 

Uncertainties in stellar parameters and e lemental abundance s 
were, as for the previous sample of 26 stars dBensby et al.l201 lh. 
calcul ated according to the method outlined in lEpstein et al.l 
J20"lol) . This method takes into account the uncertainties in the 
four observables used to find the stellar parameters, i.e. the un- 
certainty of the slope in the graph of Fei abundances versus 
lower excitation potential; the uncertainty of the slope in the 
graph of Fe 1 abundances versus reduced line strength; the un- 
certainty in the difference between Fei and Fen abundances; 
and the uncertainty in the difference between input and output 
metallicities. The method also accounts for abundance spreads 
(line-to-line scatter) as well as how the average abundances for 
each element reacts to changes in the stellar parameters (see also 
comments in lBensbv et alJl201~fl) . 

Figure|4]shows the location in the log g-T e ff plane, for all 58 
stars in the microlensed bul ge dwarf sample over-p lotted on the 



Available at |http : //obs . c arnegies cience . e du/Code/mike ^-enhanced isochrones from lPemarque et"aD d2004l) . A majority 
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Fig. 4. All 58 microlensed dwarf stars plotted on the ff-enhanced isochrones from iDemarque et~ai1 (120041) . The metallicity of the 
stars in each panel is within ±0.05 dex of the isochrone metallicity. In each plot the solid lines represent isochrones with ages of 
0.1, 5, 10, 15, and 20Gyr. Dotted lines are isochrones in steps of 1 Gyr. Error bars represent the uncertainties in T e g and \ogg as 
given in Ta ble ffl Note that the age determinations reported in Tableware based on probability distribution functions as described in 
iBensbv et all d201 ll) . 
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Fig. 5. (a) Difference between microlensing colours and the 
spectroscopic colours versus spectroscopic r e ff. (b) Difference 
between the absolute / magnitudes from microlensing tech- 
niques and from spectroscopy versus stellar mass (derived from 
spectroscopy). Error bars represent the uncertainties in the spec- 
troscopic values. The microlensing values used in the plots are 
based on the assumption that the bulge red clump has (V - I)o = 
1.06 and M, = -0.12. Stars with r eff > 5500 K are marked by 
filled circles, otherwise by open circles. 



of the stars are located either on the main-sequence turn-off or 
the subgiant branch. Three or four stars have started to ascend 
the giant branch but not far enough to have the chemical com- 
position of their atmospheres altered. The fact that most of the 
stars are either turn-off or subgiant stars, and that the uncertain- 
ties in the stellar parameters are well-constrained with usually 
small uncertainties, means that it is possible to determine rel- 
ative ages accurately. Stellar ages, masses, luminosities, abso- 
lute / magnitud es (Mi), and colours (V - 1) were estimated from 
Y 2 isochrones dDemarque et alj|2004l) b y maximising probabil- 
ity distribution functions as described in lBensbv et al] (1201 ll) . 

We have also double-checked the analysis of the previous 26 
events and revised some equivalent width measurements. This 
lead to changes in the stellar parameters and abundance ratios for 
a few stars. Also, the inclusion of NLTE corrections for Fe i lead 
to some, although small, changes. Therefore, equivalent width 
measurements and abundances for individual lines are given in 
Table[2]for all 58 stars, and stellar parameters, stellar ages, and 
elemental abundances in Tables [3] IU E) for all 58 stars, i.e., all 
tables include the revised values for the 26 stars in lBensbv et al] 
All abundances have been normalised to the Sun 
on a line -by-line basis using our own solar analysis as reference 
(see also lBensbv etal]|2011l) . 

3.2. (V - I) colours from spectroscopy and microlensing 
techniques 

The (V— 7)o colour and t he Mj magnitude can be estimated from 
microlensing techniques dYoo et al.l2004l) assuming that the red- 



Table 2. Measured equivalent widths and calculated elemental 
abundances for all 58 microlensed bulge dwarfs.' 



Element A x\ scir 1 " ' star 12 

[A] [eV] W A 6(X) ■•• W A e(X) 



' For each line we give the log cjf value, lower excitation energy (x\), 
equivalent width (W,t), absolute abundance (log £(X)). The table is only 
available in electronic form at the CDS via anonymous ftp to 
cdsarc .u-strasbg . fr (13S. 79 . 128. 5) or via 
http : //cdsweb . u- strasbg . f r/cgi -bin/qcat? J /A+A/XXX/AXX. 

Table 3. Elemental abundance ratios, errors in the abundance 
ratios, and number of lines used, for all 58 microlensed dwarf 
stars'". 



Object ... [Mg/Fe] e[Mg/Fe] N(Mg) <x(Mg) 



' For each star we give abundance ratios normalised to the Sun ([XfFe]), 
number of lines used for each element (N(X)), line-to-line scatter for each 
element (cr(X)), and the calculated uncertainty for each abundance ratio 
(e[X/X]). The table is available in electronic form at the CDS via anonymous 
ftp to cdsarc .u-strasbg . fr (138. 79. 128 . 5) or via 
http : //cdsweb . u- strasbg . f r/cgi -bin/qcat? J /A+A/XXX/AXX. 

dening towards the microlensed source is the same as towards 
the red clump in the same field, that (V - T)o and Mj of the bulge 
red clump is known, and (for Mj) that the distance to the source 
and the red clump is the same. The de-reddened magnitude and 
colour of the source can then be derived from the offsets be- 
tween the microlensing source and the red clump in the instru- 
mental colour-magnitude diagram. The microlensing values for 
Mj and (V - /) given in Table [5] are based on the assumption 
th at the bulge red clu mp has (V — T)o - 1-06 (as determine d 
in lBensbv et aDl201lh and Mj = -0.12 dNataf & Gouldll2012l) . 
For 55 of the 58 events, Fig.[5]shows comparisons between the 
spectroscopic M\ magnitudes and (V - 1) colours (see Sect. 13. U 
to th ose determined from microlensing techniques. Two events 
from iBensbv et al.ll2.Ql il and one of the new events do not have 
any estimates of (V - I)o from microlensing techniques (see 
Table [5). The first thing to notice is that there appears to be 
a slope present between the colour difference (A(V - I)) and 
the spectroscopic T e ff (see red line in Fig. |5). The question is 
whether this trend is due to uncertainties in the spectroscopic 
analysis or due to the microlensing assumptions? 

The bul ge is k nown to have patchy and irregular reddening. 
iNatafet alWM% measured reddening and differential redden- 
ing for more then 9 000 sight lines towards the bulge. Figure [6] 
shows the differential reddening along the sight lines to the in- 
dividual microlensing events in this study versus the absolute 
value of the A(V - /). While there is no one-to-one relation- 
ship between these two parameters it is evident that there are 
very few events with large colour discrepancies that at the same 
time have low differential reddening va lues (i.e., in the lower 
right corner of the Fig. ©. Note that the iNataf & Gouldl (120121) 
differential-reddening dispersions are systemically larger than 
the offsets between the colours determined from microlensing 
and spectroscopy. For example, when Nataf estimates 6E(V - 
/) ~ 0.15, the microlensing offset averages about 0.06. The 
reason for this apparent discrepancy is that to estimate the mi- 
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Fig. 6. Differential reddening versus the absolute value of 
the difference between microlensing colours and spectroscopic 
colours. The differe ntial reddening valu es for the individual sight 
lines are taken from lNataf et al.l d2012l) . Stars with r eff > 5500 K 
are marked by filled circles, otherwise by open circles. 



crolensing (V - /)o, we consider the minimum possible area 
around the sour ce needed to determine the colour of the bulge 
clump, whereas iNataf & Gouldl d2012l) use a robust algorithm 
that chooses a larger area, in order to minimise failures among 
their 9 000 automated determinations. 

In any case, Fig. [6] gives a hint that differential reddening 
could play a role in explaining the discrepancies between spec- 
troscopic and microlensing (V—I) colours. However, differential 
reddening cannot cause the (2<x statistically significant) slope in 
the A(V - 1) - T e ff diagram. The only assumptions being made 
are that the extinction is the same (on average) between the red 
clump and the source, and that the bulge clump (V - /) colour is 
known. While either of these assumptions might be violated, it 
is hard to seen why, e.g., the bulge red clump should be redder 
in fields with reddish subgiants than turnoff stars, or why there 
should be systematically more differential reddening (relative to 
the mean of that field) for fields containing reddish subgiants 
compared to turnoff stars. 

So, that turns th e focus to the spectroscopic colours. As ex- 
plained by lYi et al.l (1200 ll) . the Y 2 isochrones come with two 
dif ferent r e ff-to - colour transformation s employed: ei t her the one 
by iGreen et ail (1 19871) or the one by iLeieune et all dl998l) . As 
also pointed out bv lYi et al.l (|200l[) , while the two sets of colour 
transformations agree fairly well for main sequence dwarf stars, 
there are large differences for giant stars. The version of the 
Y 2 isochrones that we used for the mi crolensed dwa r f sam - 
ple utilises the colour transformation from lLeieune et al 1 (119981) . 
Figure [7^ shows the spectroscopic (V - /) (black filled circles) 
and the microlensing ones (open circles). Overplotted are three 
sets of r e ff-colour transformations for [ Fe/H] = and thre e 
surface gravities l ogg = 3.5, 4.0, 4.5: Leieune et al.l (Il998l) . 
IGreen et all d!987l) . and ICasagrande et all d2010l) . which have 
no dependence on log g. Figures ITJ? and c then show the dif- 
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Fig. 7. (a) Filled black circles mark the spectroscopic val- 
ues, open circles microlen s ing va lues. Dash-dotted blue line 
show the ICasagrande et alJ (20101) r e ff-co lour transformation, 
dashed green lines the Green et all dl987l) trans formation (for 
three d ifferent log g), and the solid red lines the ILeieune et al.l 
d 19981) tr ansformation (for t hree di fferen t log g). (b) Dif f erence s 
between ICasagrande et alJ d2010h and ILeieune etal] dl998l). 
c) Di fferences between IGreen et al. and ILeieune et al.l 

19981) . The markers in (b) and (c) show the differences 
that should be applied to the microle n sed d warfs if the 
ICasagrande et al.l d2010l) and IGreen et ail (ll 987b trandforma- 
tions h ad been applied instead of the one from ILeieune et aD 
(119981) . 



ferences between the Casagrande-Lejeune and Green-Lejeune 
transformations, respectively. What is interesting here is that the 
Lejeune and Casagrande transformations are nearly identical for 
T e ff > 5500 K, and that there is an offset between Green and 
Lejeune that varies with log g but which is nearly constant for 
Jeff > 5500 K. However, for r e jf <5500 the three sets of transfor- 
mations differ, and increasingly more for lower T e g . The points 
in Figs.|7p and c show the corrections that would have to be ap- 
plied to the microlensed dwarf stars if we were to change from 
the Lejeune transformation to any of the other two. Generally, 
stars with r e ff > 5500 would be more or less untouched while the 
cooler stars would get higher {V—I) colours. This is a frustrating 
finding as it means that the slope seen in Fig. [5] would become 
steeper with the Casagrande and Green r e ff-colour transforma- 
tions ! 

However, from Fig. [7] it is apparent that the agreement be- 
twe en the three calibrat ions is good for hotter stars. Even though 
the IGreen etal] dl987l) transformations will show offsets de- 
pending on log g, it will be a constannt offset and no slope will be 
introduced. As further investigations into the field of r e ff-colour 
transformations is beyond the scope of the current paper, we will, 
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for now, restrict ourselves to stars with T e ff > 5500 K when com- 
paring spectroscopic and microlensing colours. These stars are 
marked by filled circles in Fig. [5] The difference in the (V - I) 
colour is now 0.00 + 0.07 mag (while if we were to include all 
stars, i.e., also those with < 5500 K, the difference would be 
+0.03 + 0.08 mag). Our conclusion is that the colour of the bulge 
red clump should remain at (V - 7)o = 1 .06. 

3.3. Comparisons to the Balmer Ha line 

To further check the effective temperatures of the microlensed 
dwarf stars we have calculated synthetic sp ectra with the SME 
(Spec troscopy Made Easy, v. 201 1-12-05, [Valenti & Piskunovl 
1 19961) for the 32 new microlensing events. Figure [8] shows 
comparisons of the synthetic Ha line profiles to the observed 
Ha line profiles the spectroscopic temperatures (blue lines). 
There is generally very good agreement between synthetic and 
observed spectra for the spectroscopic effective temperatures. 
Hence, we believe that the effective temperatures we have de- 
termined should be good. 

For the 55 stars that have (V - I)q colours from microlens- 
ing techniques we calc ulate effective temperatures using the 
ICasagrande et al.l d2010l) r e ff-colour calibration. These temper- 
atures are on average 105 K higher than the spectroscopic ones. 
If we restrict ourselves to stars with (spectroscopic) T e ff > 5500 
the spectroscopic temperatures are on average higher, but only 
by 10K. The dispersion in both cases is around 230 K. 

Figure [8] also shows synthetic spectra based on tempera- 
tures from the microlensing (V - 7)o colours (red dashed lines). 
It is clear that the microlensing temperatures produce spectra 
that do not match for a few cases (e.g., MOA-201 1-BLG-090S 
and MOA-2012-BLG-410S). The Ha wing profile for the star 
with the largest colour discrepancy between spectroscopic and 
microlen sing temperature, MO A-2009-BLG-259S, is shown in 
Fig. 3 of lBensbv etaTI d20lTl) . From that figure it is clear that 
the spectroscopic temperature is the better match. However, this 
star has r e jf = 4915 K and loggr = 3.3, which is in a region 
where the r e ff-colour calibrations appear to be very uncertain 
(see Fig. [7J, meaning that the colour- T e $ transformation could 
be substantially affected. 

3.4. A few special cases 

The combination of low metallicity and low S /N of the ob- 
served spectra of OGLE-201 1-BLG-0969S and OGLE-20 12- 
BLG-1156S meant that very few equivalent widths could be 
measured for these stars. Hence, the effective temperatures and 
the microturbulence parameters for these two stars could not 
be constrained in the standard way, through excitation and ion- 
isation balance, as described above. As an alternative, the ef- 
fective temperatures were constrained by fitting the wing pro- 
files of the Balmer Ha line at 6563 A for these two stars (as 
described in Sect. 13.2b . This gives an effective temperature of 
around 6 1 50 K for OGLE-20 1 1 -BLG-0969S and around 6200 K 
for OGLE-20 1 2-BLG- 1 156S (see Fig.|9]l. This in fair agreement 
with the temperatures that can be estimated from microlens- 
ing techniques for these two stars, 6283 K and 6012 K, respec- 
tively (see Sect. 13.21 1. The surface gravities were set based on 
microlensing techniques to \ogg = 4.14 for OGLE-201 1-BLG- 
0969S and logg = 4.25 for OGLE-20 1 2-BLG- 1 156S. First, 
(V - 7)o and Io are determined from the offset to the bulge red 
clump. (V - I)p is then transferre d to (V - K)o using the cal- 
ibration by iBessell & Bretl dl988l) . Combining (V - K) with 
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Fig. 9. Fitting of the Ha wing profiles to find the effective tem- 
peratures for OGLE-201 1-BLG-0969S and OGLE-20 12-BLG- 
1 156S. The insets show the summed absolute values of the resid- 
uals that was minimised to find the best fitting temperature. 

Vq - ( V-I)o+Iq give slar using the calibration bv lKervella et all 
d200l . We then assume that the mi crolensing source i s at same 
distance as the red clump given by iNataf et all d2012l) . The ra- 
dius of the star is then r star = D source ■ 9 slar . Finally, we as- 
sume that the stellar mass is M slar = 1 M , so the gravity is 
given by g — GM slar lr 2 slar We estimate the uncertainty in the 
distance to be 15% (0.3 mag), and the mass uncertainty to be 
about 10%. Hence, the dominant error is in the distance (because 
it enters as the square), imply that the error in log g using this 
method is about 0. 1 dex. The microturbulence parameter was set 
to £ t = L5 kms~' for both stars using the empirical calibration 
by M. Bergemann (private communication^ 

Following the standard procedure to determine the stellar pa- 
rameters (see above) for MOA-201 1-BLG-104S we get a mi- 
croturbulence value of £ t = 2.7 kms . This value seems too 
high for these types of stars. According to the empirical cali- 
bration by M. Bergemann, stars with similar stellar parameters 
{T eS = 6106 K, log g = 4.4, and [Fe/H] = -0.88) have mi- 
croturbulence values around 1.5 kms -1 . Also, the discrepancy 
between the spectroscopic temperature and the temperature de- 
duced from microlensing techniques is unusually large for this 
star (~ 500 K). Hence, we believe that due to the quite low metal- 
licity of this star, combined with a relatively low signal-to-noise 
ratio of the observed spectrum (S IN « 25), the stellar parame- 
ters have been poorly constrained from the Fe i balance plots. 

The discrepancies indicate that it is likely that the temper- 
ature should be lower than the 6100K derived from excitation 
balance. However, based on the Balmer Ha line profile it should 
probably not be as low as the temperature deduced from the mi- 
crolensing techniques (5644 K, see Table [5] and Fig. [8}. Hence, 

2 The empirical microturbulence relation by Maria B ergemann was 
deriv ed for the Gaia-ESO public spectroscopic survey iGilmore et al .1 
120121) and will be published in one of her first science papers based on 
Gaia-ESO data. 
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6555 6560 6565 6570 
Wavelength [A] 



6555 6560 6565 6570 
Wavelength [A] 



6555 6560 6565 6570 
Wavelength [A] 



Fig. 8. Illustration to see how well the spectroscopic temperatures (full blue lines) and temperatures from microlensing techniques 
(dashed red lines) reproduce the observed wing profiles of the Ha line at 6563 A for the 32 new microlensed bulge dwarfs. Note 
that no estimate of the temperature from microlensing techniques is available for OGLE-201 1-BLG-1410S. The stars are sorted by 
metallicity. 
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Table 4. Stellar parameters and radial velocities for all 58 stars (sorted by metallicity). 1 



Object 


T e g 


logs 


ft 


[Fe/H] 


Ate I, Fen 






[K] 


[cgs] 


[kms- 1 ] 






[km/s] 


OGLE-2012-BLG-1156S 


6200 ± 200 


4.25 ± 0.30 


1.50 ±0.30 


-1.89 ±0.25 


3, 


72.6 


OGLE-2011-BLG-0969S 


6150 ± 200 


4.14 ±0.30 


1.50 ±0.30 


-1.57 ±0.25 


11, 2 


-53.0 


MOA-2010-BLG-285S 


6057 ± 134 


4.20 ±0.19 


1.88 ±0.42 


-1.21 ±0.10 


53, 9 


46.0 


MOA-2010-BLG-078S 


5205 ± 189 


3.60 ± 0.26 


1.31 ±0.32 


-1.00 ±0.34 


53, 4 


52.3 


MOA-2011-BLG-104S 


5900 ± 200 


4.15 ±0.30 


1.30 ±0.30 


-0.85 ± 0.25 


36, 6 


197.3 


OGLE-20 1 2-BLG-0270S 


5914 ± 145 


4.30 ± 0.22 


1.33 ±0.22 


-0.84 ±0.13 


45, 6 


-128.7 


MOA-2009-BLG-493S 


5420 ± 119 


4.40 ± 0.22 


0.90 ± 0.32 


-0.74 ±0.15 


80, 5 


-14.5 


MOA-2012-BLG-187S 


5895 ± 94 


4.20 ±0.14 


1.60 ±0.18 


-0.74 ± 0.08 


72, 8 


-40.3 


OGLE-2009-BLG-076S 


5854 ± 108 


4.30 ±0.15 


1.63 ±0.22 


-0.72 ±0.11 


57, 7 


128.7 


MOA-2009-BLG-133S 


5529 ± 73 


4.30 ±0.14 


1.17 ± 0.18 


-0.69 ± 0.07 


66, 6 


91.6 


OGLE-2012-BLG-0563S 


5907 ± 89 


4.40 ±0.10 


1.27 ±0.14 


-0.66 ± 0.07 


73, 8 


-66.2 


OGLE-20 12-BLG-1279S 


5796 ± 63 


4.40 ±0.11 


0.99 ±0.12 


-0.62 ± 0.06 


152, 18 


141.1 


MOA-2010-BLG-167S 


5406 ± 49 


3.90 ± 0.09 


1.15 ±0.09 


-0.60 ± 0.05 


109, 14 


-79.4 


MOA-2012-BLG-532S 


5626 ± 207 


3.90 ±0.37 


0.88 ± 0.25 


-0.55 ±0.21 


42, 2 


27.9 


MOA-2009-BLG-475S 


5836 ± 189 


4.40 ± 0.27 


1.35 ±0.37 


-0.52 ± 0.20 


53, 4 


137.8 


MACH- 1 999-BLG-022S 


5632 ± 110 


4.10±0.18 


0.33 ±0.41 


-0.49 ±0.17 


97, 10 


37.6 


OGLE-2012-BLG-1217S 


5795 ± 73 


4.30 ±0.13 


1.28 ±0.15 


-0.41 ±0.07 


73, 9 


124.4 


MOA-2010-BLG-446S 


6308 ± 111 


4.50 ±0.14 


1.71 ±0.18 


-0.40 ± 0.08 


66, 8 


56.5 


MOA-2010-BLG-049S 


5694 ± 61 


4.10±0.12 


1.02 ±0.12 


-0.40 ± 0.07 


96, 14 


-116.7 


OGLE-2008-BLG-209S 


5248 ± 77 


3.80 ±0.14 


1.08 ±0.10 


-0.31 ±0.09 


147, 18 


-173.6 


MOA-2011-BLG-090S 


5367 ± 49 


4.10 ±0.09 


0.87 ± 0.09 


-0.26 ± 0.05 


88, 8 


48.2 


MOA-2012-BLG-391S 


5505 ± 76 


3.90 ±0.14 


1.12 =t 0.16 


-0.24 ± 0.07 


67, 7 


-65.0 


OGLE-20 12-BLG-1526S 


5200 ± 62 


3.85 ±0.12 


0.94 ±0.12 


-0.24 ± 0.06 


76, 6 


-87.1 


MOA-2009-BLG-489S 


5543 ±61 


4.10 =t 0.1 1 


0.62 ±0.10 


-0.21 ±0.07 


106, 12 


96.5 


MOA-2011-BLG-174S 


6172 ± 111 


4.40 ±0.16 


1.16 ± 0.14 


-0.18 ±0.09 


57, 6 


-24.0 


MOA-2012-BLG-202S 


4862 ± 93 


3.90 ±0.18 


0.90 ±0.17 


-0.15 ±0.13 


60, 4 


41.2 


OGLE-2012-BLG-1534S 


5920 ± 52 


4.05 ± 0.08 


1.37 ±0.10 


-0.15 ±0.04 


94, 12 


206.4 


MOA-2012-BLG-410S 


5509 ± 80 


3.90 ±0.13 


1.09 ±0.14 


-0.14 ±0.07 


60, 7 


22.3 


OGLE-20 12-BLG-0617S 


4924 ±71 


3.70 ±0.16 


1.06 ±0.11 


-0.14 ±0.09 


85, 8 


-68.0 


OGLE-20 12-BLG-0816S 


5395 ± 131 


3.50 ± 0.20 


0.65 ± 0.20 


-0.10 ±0.14 


37, 5 


-74.0 


OGLE-2012-BLG-0211S 


5573 ±75 


4.00 ±0.12 


0.67 ±0.13 


-0.06 ± 0.08 


79, 10 


-17.7 


MOA-2011-BLG-234S 


5778 ± 88 


4.40 ±0.13 


0.93 ± 0.08 


-0.02 ± 0.08 


111, 15 


44.0 


OGLE-2011-BLG-1105S 


5666 ± 113 


4.50 ±0.18 


0.75 ± 0.22 


+0.00 ±0.17 


37, 5 


222.7 


MOA-2010-BLG-523S 


5122 ±79 


3.60 ±0.15 


1.68 ±0.20 


+0.06 ±0.14 


57, 10 


97.3 


OGLE-20 12-BLG-1274S 


5733 ±51 


4.10 ±0.07 


1.23 ±0.07 


+0.07 ± 0.04 


94, 12 


-25.0 


OGLE-2012-BLG-0521S 


5013 ±84 


3.70 ±0.14 


1.02 ±0.11 


+0.09 ±0.15 


75, 10 


-68.8 


MOA-2009-BLG-174S 


5600 ± 80 


4.40 ±0.10 


1.27 ±0.10 


+0.11 ±0.08 


117, 15 


-21.4 


MOA-2011-BLG-034S 


5440 ± 91 


4.10±0.12 


0.88 ±0.11 


+0.11 ±0.15 


77, 10 


127.0 


MOA-2009-BLG-456S 


5662 ± 89 


4.20 ±0.11 


0.86 ±0.10 


+0.13 ±0.10 


88, 12 


-164.6 


MOA-2012-BLG-291S 


5156 ± 107 


4.10 ±0.15 


0.79 ±0.13 


+0.16 ±0.24 


64, 9 


60.3 


OGLE-2011-BLG-1410S 


4831 ± 108 


3.30 ± 0.24 


0.62 ±0.15 


+0.22 ± 0.37 


37, 4 


-75.2 


MOA-2011-BLG-191S 


5382 ± 92 


3.80 ±0.13 


0.57 ±0.11 


+0.26 ±0.14 


61, 11 


134.0 


MOA-2011-BLG-445S 


4870 ± 102 


3.50 ±0.19 


0.61 ±0.12 


+0.26 ±0.31 


57, 8 


72.6 


OGLE-2007-BLG-514S 


5644 ± 111 


4.10 ±0.21 


1.47 ±0.21 


+0.29 ± 0.23 


43, 7 


158.8 


OGLE-2011-BLG-0950S 


6130 ± 121 


4.20 ±0.15 


1.23 ±0.15 


+0.33 ±0.10 


58, 7 


91.5 


MOA-2009-BLG-259S 


4915 ± 104 


3.30 ±0.18 


1.01 ±0.10 


+0.34 ±0.19 


64, 8 


81.5 


MOA-2008-BLG-311S 


5947 ±81 


4.50 ±0.10 


1.17 ±0.09 


+0.35 ± 0.08 


111, 14 


-34.1 


OGLE-2011-BLG-1072S 


5515 ±89 


3.90 ±0.13 


1.11 ±0.10 


+0.36 ±0.12 


69, 13 


-62.2 


MOA-2011-BLG-058S 


5256 ± 100 


4.00 ±0.15 


0.71 ±0.12 


+0.37 ± 0.25 


62, 10 


-139.7 


MOA-2008-BLG-3 10S 


5675 ± 91 


4.20 ±0.11 


1.03 ± 0.08 


+0.41 ±0.11 


119, 18 


77.5 


MOA-2012-BLG-022S 


5827 ± 115 


4.30 ±0.13 


0.89 ±0.12 


+0.42 ±0.10 


53, 9 


-81.3 


OGLE-2007-BLG-349S 


5237 ± 119 


4.20 ±0.15 


0.77 ±0.12 


+0.42 ± 0.26 


103, 17 


113.0 


MOA-2006-BLG-099S 


5747 ± 147 


4.50 ±0.16 


0.84 ±0.14 


+0.44 ±0.21 


117, 17 


99.0 


OGLE-2006-BLG-265S 


5462 ± 101 


4.20 ±0.13 


1.14 ± 0.11 


+0.46 ±0.19 


92, 14 


-154.0 


OGLE-20 1 2-BLG-0026S 


4815 ± 145 


3.40 ± 0.28 


0.62 ±0.14 


+0.50 ± 0.44 


42, 7 


132.2 


MOA-2010-BLG-311S 


5442 ± 86 


3.80 ±0.15 


1.16 ± 0.11 


+0.51 ±0.19 


52, 11 


44.4 


MOA-2011-BLG-278S 


5307 ± 159 


4.00 ± 0.25 


0.54 ±0.18 


+0.52 ± 0.39 


48, 7 


229.6 


MOA-2010-BLG-037S 


5732 ± 109 


3.90 ± 0.20 


1.52 ±0.19 


+0.55 ± 0.20 


56, 12 


-8.4 



Notes. ' Given for each star is: effective temperature (T e fi ); surface gravity (log g); microturbulence parameter (if,); [Fe/H]; number of Fe i and Fe n lines used in 
analysis; stellar age; radial velicity (t) r ). 

for MOA-2011-BLG-104S, we set the temperature to 5900 K. -0.85). The age of the star increased from 8.7 Gyr to 10.2 

Following the empirical calibration by M. Bergemann we set (which is more in line with ages seen for the other metal-] 

the microturbulence to 1.3 kms~'. To get ionisation balance be- stars). 

tween Fe i and Fe n abundances the surface gravity needs to be 

set to loggr = 4.15. The effect on the metallicity for MOA- 

201 1-BLG-104S, when changing the temperature from 6100K 

to 5900K is only -0.03 dex (down from [Fe/H] = -0.88 to 
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Table 5. Stellar ages, masses, and colours as determined from spectroscopy and microlensing techniques for all 58 stars. 1 



WDJcCI 


7 

' efl 


A 1 


VOgL 


Age 


1 rr 
— LU 


+ 1 (T 


\y-')o 


AA 


1 ,T 


+ 1 (T 


Ml 




(' 1 10,11 


...ij 

'ell 




[K] 


M.-.j 


Lq 


[Gyr] 


[Gyr] 


IGyrJ 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[K] 


OGLE-20 1 2-BLG- 1156 


6200 


0.77 


0.05 


13.1 


6.4 


13.4 


0.58 


4.79 


3.43 


5.19 


4.21 


3.56 


0.63 


6058 


OGLE-201 l-BLG-0969 


6150 


0.79 


0.02 


11.6 


5.9 


13.2 


0.59 


4.86 


3.85 


5.10 


4.27 


2.94 


0.58 


6296 


MOA-20 1 0-BLG-285 


6057 


0.78 


0.09 


11.6 


8.2 


13.4 


0.62 


4.67 


3.86 


4.99 


4.05 


4.77 


0.56 


6384 


MOA-2010-BLG-078 


5205 


0.93 


0.77 


14.6 


4.9 


13.1 


0.82 


2.99 


2.38 


3.57 


2.17 


1.79 


0.83 


5317 


MOA-20 11-BLG- 104 


5900 


0.81 


0.00 


10.2 


5.5 


12.9 


0.65 


4.92 


3.48 


5.20 


4.27 


3.50 


0.73 


5644 


OGLE-20 1 2-BLG-0270 


5914 


0.84 


-0.07 


11.6 


6.1 


13.0 


0.65 


5.07 


4.13 


5.25 


4.42 


3.07 


0.62 


6081 


MOA-2009-BLG-493 


5420 


0.72 


-0.40 


13.1 


4.7 


13.3 


0.79 


5.98 


5.59 


6.20 


5.19 


4.74 


0.68 


5829 


MOA-20 1 2-BLG- 187 


5895 


0.84 


0.19 


11.6 


9.2 


13.4 


0.67 


4.40 


4.02 


4.87 


3.73 


3.71 


0.71 


5715 


OGLE-2009-BLG-076 


5854 


0.83 


-0.03 


11.6 


7.6 


13.4 


0.68 


4.96 


4.32 


5.18 


4.28 


4.17 


0.68 


5828 


MOA-2009-BLG-133 


5529 


0.75 


-0.31 


14.6 


7.4 


13.9 


0.76 


5.71 


5.36 


5.81 


4.95 


4.17 


0.69 


5790 


OGLE-20 12-BLG-0563 


5907 


0.85 


-0.04 


10.2 


5.7 


12.0 


0.66 


4.89 


4.63 


5.13 


4.23 


4.36 


0.58 


6257 


OGLE-20 1 2-BLG- 1 279 


5796 


0.83 


-0.07 


11.6 


6.6 


13.1 


0.69 


5.06 


4.76 


5.26 


4.37 


4.11 


0.61 


6117 


MOA-20 10-BLG- 167 


5406 


0.87 


0.37 


13.5 


9.0 


13.8 


0.80 


4.04 


3.65 


4.07 


3.24 


2.71 


0.70 


5751 


MOA-20 12-BLG-532 


5626 


0.92 


0.37 


11.6 


4.2 


12.7 


0.75 


3.99 


3.00 


5.01 


3.24 


4.15 


0.67 


5864 


MOA-2009-BLG-475 


5836 


0.87 


-0.05 


8.7 


4.4 


12.3 


0.69 


5.04 


4.19 


5.37 


4.35 


4.23 


0.60 


6157 


MACH- 1 999-BLG-022 


5632 


0.85 


0.29 


13.1 


8.0 


13.6 


0.74 


4.17 


3.66 


4.85 


3.43 








OGLE-20 1 2-BLG- 1217 


5795 


0.89 


0.00 


11.6 


7.2 


12.4 


0.69 


4.76 


4.35 


5.03 


4.07 


3.67 


0.71 


5712 


MOA-20 10-BLG-049 


5694 


0.87 


0.25 


12.4 


9.4 


13.6 


0.73 


4.25 


3.82 


4.55 


3.52 


2.67 


0.67 


5861 


MOA-20 10-BLG-446 


6308 


1.06 


0.20 


2.6 


1.3 


4.4 


0.57 


4.32 


4.00 


4.51 


3.75 








OGLE-2008-BLG-209 


5248 


0.97 


0.40 


8.6 


5.3 


12.5 


0.84 


3.92 


3.34 


4.09 


3.08 


2.50 


0.74 


5608 


MOA-20 ll-BLG-090 


5367 


0.91 


0.26 


14.3 


11.1 


14.4 


0.80 


4.30 


4.02 


4.30 


3.50 


3.77 


0.67 


5857 


MOA-20 12-BLG-391 


5505 


1.00 


0.40 


8.1 


5.3 


12.2 


0.77 


3.89 


3.33 


4.15 


3.12 


3.92 


0.81 


5390 


OGLE-20 1 2-BLG- 1 526 


5200 


0.93 


0.37 


9.6 


6.4 


13.0 


0.84 


3.99 


3.54 


4.17 


3.15 


2.99 


0.69 


5782 


MOA-2009-BLG-489 


5543 


0.92 


0.26 


14.2 


9.8 


14.0 


0.76 


4.25 


3.91 


4.55 


3.49 


3.37 


0.86 


5251 


MOA-201 l-BLG-174 


6172 


1.08 


0.18 


3.1 


1.7 


5.2 


0.60 


4.37 


3.85 


4.50 


3.77 


4.27 


0.66 


5893 


MOA-20 1 2-BLG-202 


4862 


0.97 


0.41 


14.7 


8.6 


13.9 


0.94 


3.96 


3.56 


4.17 


3.02 


5.09 


0.90 


5151 


OGLE-20 1 2-BLG- 1 534 


5920 


1.08 


0.44 


6.8 


5.8 


7.7 


0.65 


3.73 


3.44 


3.94 


3.08 


2.61 


0.64 


5971 


MOA-20 12-BLG-4 10 


5509 


1.01 


0.46 


6.8 


5.2 


11.7 


0.77 


3.73 


3.33 


4.13 


2.96 


5.32 


0.94 


5056 


OGLE-20 1 2-BLG-06 1 7 


4924 


0.95 


0.48 


13.3 


7.0 


13.4 


0.94 


3.79 


3.30 


4.02 


2.85 


3.42 


1.03 


4867 


OGLE-20 12-BLG-08 16 


5395 


1.29 


0.81 


1.6 


1.7 


7.0 


0.80 


3.30 


1.98 


3.63 


2.50 


3.20 


0.72 


5673 


OGLE-20 1 2-BLG-02 1 1 


5573 


1.00 


0.35 


8.6 


6.5 


11.9 


0.77 


3.98 


3.52 


4.31 


3.21 


2.99 


0.70 


5742 


MOA-201 l-BLG-234 


5778 


0.99 


0.00 


6.4 


3.3 


8.8 


0.69 


4.76 


4.44 


4.98 


4.07 


4.98 


0.74 


5606 


OGLE-201 1-BLG- 1105 


5666 


0.95 


-0.07 


6.4 


3.1 


10.4 


0.74 


5.05 


4.62 


5.26 


4.31 


4.03 


0.78 


5480 


MOA-20 1 0-BLG-523 


5122 


1.14 


0.58 


3.3 


2.8 


8.8 


0.90 


3.53 


2.80 


3.90 


2.63 




0.75 


5573 


OGLE-20 1 2-BLG- 1 274 


5733 


1.05 


0.32 


8.3 


7.2 


9.2 


0.71 


4.03 


3.75 


4.21 


3.32 


3.39 


0.68 


5811 


OGLE-20 12-BLG-0521 


5013 


1.03 


0.47 


5.7 


4.6 


11.8 


0.91 


3.82 


3.30 


4.13 


2.91 


2.98 


0.85 


5283 


MOA-2009-BLG- 1 74 


5600 


0.96 


-0.05 


8.7 


4.2 


10.4 


0.76 


5.01 


4.70 


5.18 


4.25 


3.50 


0.69 


5774 


MOA-201 l-BLG-034 


5440 


0.94 


0.23 


12.4 


8.7 


13.6 


0.79 


4.31 


3.92 


4.67 


3.52 




0.65 


5922 


MOA-2009-BLG-456 


5662 


1.01 


0.14 


8.7 


6.7 


10.9 


0.74 


4.39 


4.06 


4.74 


3.65 


2.74 


0.67 


5846 


MOA-20 12-BLG-291 


5156 


0.96 


0.27 


14.6 


8.2 


13.8 


0.87 


4.29 


3.94 


4.67 


3.42 


4.02 


0.91 


5135 


OGLE-201 l-BLG-1410 


4831 


1.02 


0.86 


6.5 


3.4 


12.1 


1.00 


2.91 


2.20 


3.51 


1.91 


2.13 






MOA-201 l-BLG-191 


5382 


1.18 


0.50 


3.9 


3.7 


9.2 


0.82 


3.65 


3.08 


4.00 


2.83 


2.50 


0.67 


5841 


MOA-201 l-BLG-445 


4870 


0.98 


0.68 


8.2 


4.1 


12.1 


0.98 


3.32 


2.78 


3.86 


2.34 


3.20 


0.92 


5114 


OGLE-2007-BLG-5 14 


5644 


1.01 


0.22 


6.4 


4.6 


10.3 


0.75 


4.33 


3.58 


4.79 


3.58 


4.34 


0.71 


5700 


OGLE-201 l-BLG-0950 


6130 


1.27 


0.34 


2.9 


1.6 


3.6 


0.60 


3.95 


3.26 


4.14 


3.35 


3.41 


0.63 


5990 


MOA-2009-BLG-259 


4915 


1.17 


0.88 


3.0 


2.0 


9.0 


1.00 


2.84 


2.06 


3.36 


1.84 


2.88 


0.80 


5421 


MOA-2008-BLG-311 


5947 


1.17 


0.13 


1.7 


0.9 


3.0 


0.66 


4.49 


4.25 


4.60 


3.83 


3.91 


0.67 


5838 


OGLE-201 l-BLG-1072 


5515 


1.18 


0.47 


5.5 


4.4 


8.8 


0.79 


3.70 


3.19 


4.08 


2.91 


2.66 


0.79 


5450 


MOA-201 l-BLG-058 


5256 


0.97 


0.26 


10.8 


6.5 


13.1 


0.85 


4.30 


3.74 


4.52 


3.45 


3.29 


0.80 


5421 


MOA-2008-BLG-310 


5675 


1.12 


0.27 


4.9 


4.4 


8.0 


0.75 


4.17 


3.86 


4.58 


3.42 


3.44 


0.70 


5731 


MOA-20 12-BLG-022 


5827 


1.15 


0.17 


3.5 


2.1 


5.5 


0.70 


4.52 


3.93 


4.65 


3.82 


3.70 


0.82 


5366 


OGLE-2007-BLG-349 


5237 


0.94 


-0.13 


13.1 


8.1 


13.6 


0.85 


4.37 


4.28 


5.47 


3.52 


4.14 


0.79 


5450 


MOA-2006-BLG-099 


5747 


1.06 


0.02 


3.5 


1.8 


7.1 


0.73 


4.82 


4.39 


5.06 


4.09 


3.79 


0.75 


5568 


OGLE-2006-BLG-265 


5462 


1.00 


0.10 


8.7 


6.3 


11.8 


0.81 


4.64 


4.19 


5.00 


3.83 


3.57 


0.69 


5764 


OGLE-20 12-BLG-0026 


4815 


1.02 


0.74 


6.5 


3.7 


12.4 


1.00 


3.15 


2.36 


3.80 


2.15 


2.68 


1.04 


4878 


MOA-20 10-BLG-3 11 


5442 


1.20 


0.50 


4.4 


3.3 


8.0 


0.81 


3.63 


2.89 


3.96 


2.82 


3.14 


0.73 


5630 


MOA-201 l-BLG-278 


5307 


1.01 


0.34 


10.2 


5.2 


12.7 


0.82 


4.09 


3.59 


5.06 


3.27 


4.49 


0.86 


5263 


MOA-20 10-BLG-037 


5732 


1.29 


0.56 


4.2 


3.1 


6.3 


0.72 


3.44 


2.85 


4.01 


2.72 


3.54 


0.78 


5478 



Notes. ' Column 2 effective temperature; col. 3 stellar mass; col. 4 luminosity, col. 5 stellar age; cols. 6 and 7 the 1-sigma lower and upper age limits; col. 8 
"spectroscopic" colour, col. 9 "spectroscopic" My magnitude, cols. 10 and 1 1 the 1-sigma lower and upper limits on My; col. 12 "spectroscopic" M\ magnitude, 
col. 13 Mi ^ from microlensing techniques (with the bulge red clump at -0.12); col. 14 (V — /) p u from mic rolensing techniques (with the bulge red clump at 1.06); 
col. 15 inferred effective temperatures from microlensing (V - /) o colour using the Casagrande et al. (2010) (V* - /) - T s g calibration. 



4. Results 

4.1. Positions on the sky 

Figure [10] shows the positions on the sky for the 58 microlensed 
dwarf stars. Except for the first three detections at positive lat- 



itudes, essentially all stars are located between 2° to 5° below 
the the Galactic plane and mainly at positive longitudes. This 
distribution is a reflection of the locations of the fields surveyed 
by MOA and OGLE. At the distance of the bulge, 8.2 kpc (e.g., 
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Fig. 10. Positions on the sky for the microlensed dwarf sample. 
The bulge contour lines based on observations with t he COBE 
satellite are shown as solid lines dWeiland et al.lll994l) . The po- 
sitions of Baade's window (BW), the Galactic centre (G C), and 
three of the ARGOS survey fields (see lNess et alJl2012h are de- 
lineated. The dotted lines are concentric circles in steps of 2°. 



iNataf et al.ll20"l2h . 2° to 5° corresponds to a vertical distance of 
150 to 750 pc, respectively, below the plane. 

4.2. Radial velocities - bulge membership 

Radial velocities can help constra ining dynamical mo dels of the 
bulge. The BRAVA survey ( see. IKunder et al.1 12Q12L for latest 
release) presents radial velocities for ~ 10000 red giant stars 
covering -10 < I < 10 and -8 < b < -4. They find that the 
rotation curves at / = -4, -6 a nd -8 almost perfe ctly follow the 
cylindrical rotation models by IShen et aT | dioiOt) . The velocity 
dispersion varies from around 100 kms -1 at (I, b) = (0, -4) to 
around 80 kms" 1 at (I, b) = (±10, -4). 

We correct our measured heliocentric radial velocities for the 
microlensed bulge dwarfs to the Galactic centre using the rela- 
tiotfl vgc = v r nc + 232 sin(Z) cos(b) + 9 cos(/) cos(fr) + 7 sin(fo) 
which takes the motion of the relative to the local standard of rest 
(LSR), and the motion of the LSR relative to the Galactic cen- 
tre, into account. Figure QT| then shows the galactocentric radial 
velocities versus I and b for the microlensed dwarf sample (grey 
circles) as well as the results from the BRAVA survey in differ- 
ent / and b bins (as indicated in the figure; red, cyan, and green). 
The black circles mark the mean velocity and the rms disper- 
sion around the mean the microlensed dwarf stars, binned into 
3° wide bins. The average velocities as well as the dispersions 
for the microlensed dwarf stars follow nicely the results from 
the BRAVA survey, especially versus I. Given that the BRAVA 
sample is representative of the bulge population, this shows that 
the radial velocities and the velocity dispersion seen in the mi- 
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Fig. 11. Galactocentric velocity versus Galactic longitude and 
latitude for the microlensed dwarf sample (grey circles). The 
black circles represent the mean radial velocity and the l<x rms 
dispersion for the microlensed bulge dwarfs, bi nned into 3° wide 
bins. The results f rom the BRAVA survey dRich et all 120071: 
IKunder et al.l 120121) have been over-plotted, and the different 
colours (red, cyan, and green) represent the values at different 
/ and b as given in the plots. 



crolensed dwarf sample are at the expected levels if they belong 
to a bulge population. 



4.3. Metallicity distribution 

FigurefTZh shows the metallicity distribution for the mi crolensed 
dwarf sample. The finding from iBensbv et al.l (1201 ll) . that the 
MDF is bi-modal, with a paucity of stars around solar metal- 
licity, is weakened as we now have several stars with metallici- 
ties within +0.1 dex of the Sun. One o f those stars, MOA-2 010- 
BLG-523S, was already published in IBensbv etail (1201 lh but 
has now been revised 0.08 dex lower metallicity following the 
current analysis (see Sect. 13. U . 

A major issue in IBensbv et al.1 (1201 ll) was the apparent dis- 
crepancy between the microlensed dwarf star MDF, which was 
clearly bi-modal, an d the red giant MDF in Baade's window by 
IZoccalietal1(l2008h. which then did not appear bi-modal. Since 
th en, thelZoccali et al.l d2008) RGB sample has been re-analysed 
bv lHill et al.1 d201 lh with the result that the RGB sample shrank 
from 213 stars to 166 stars and that the MDF now appears bi- 
modal (see Fig. [12b). With the current microlensed dwarf sam- 
ple, and the re-analysed RGB sample, the two distributions are 
quite similar. A two-sided KS test yields a />value of 0.17 (see 
Fig. [T3l , i.e., the hypothesis that the two samples have been 
drawn from the same underlying population cannot be rejected. 

et al.l (1201 ll) also presented a sample of 220 red clump 
giants in Baade's window. The MDF of these red clump stars 
(Fig. [12b ) is similar to the RGB sample, and a two-sided KS test 
between the red clump sample and the microlensed dwarf sam- 
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Fig. 12. (a) Metallicity distribution for the microlensed dwarf 
sample (white dashed line shows the distributi on of the 26 mi- 
crolensed dwarf stars from lBensbv et al.ll2.Ql 1 ); (b) 1 66 red gi- 

ic) 220 



ant stars in Baade's window from iHill et al 



red clump stars in Baade's window from IHill 



et all (l2i 



fcoill) : (d) 

1813 red giant stars with from the ARGOS survey fields at 
(l,b) = (0,-5), (5,-5), (-5,-5) from Ness e t al. (submitted) 
(e) 36 3 red giants at (l,b) = (0,-10) from lUttenthaler et al] 
(12011 : The curved lines in (a)-(d) represent generalised his- 
tograms. Dotted vertical lines mark the peaks claimed by Ness 
et al. in (d). 



pie yields a /7-value of 0.06 (see Fig.QjJi. Hence, the conclusions 
are the same as for the RGB sample. 

Furthermore, ARGOS (Abundances and Radial velocity 
Galactic Origins Survey) observed 28,000 stars, mostly red 
clump giants, in 28 fields of the Galactic bulge. In Fig. fT2h we 
show the MDF for the 1813 red clump giants that are located 
in the three ARGOS fields at {l,b) = (0, -5), (5, -5), (-5, -5) 
(the fields are delineated out in Fig. HoV The internal ran- 
dom uncertainty in [Fe/H] for individual stars in ARGOS is 
0.13dex (Ness et al. submitted). By fitting Gaussians to the 
MDF, Ness et al. claim detection of three components at 
[Fe/H] = +0.14, -0.23, -0.60, and possibly another two at 
[Fe/H] = -1.24, -1.7. The peaks are indicated by the vertical 
dotted lines in Fig. [12] It is interesting to see that the three metal- 
rich ARGOS peaks more or less coincide with the "bumps" in 
the generalised dwarf star MDF in Fig. [T2k . There might be 
a slight offset, but the three peaks are clearly within 0. 1 dex 
of the "bumps". It is also interesting to note that the ARGOS 
MDF does not show the relatively large fraction of metal-rich 
stars ([Fe/H] > +0.3) seen in the microlensed dwarf sample 
and the RGB and RC samples in Baade's window. Since the 
ARGOS fields are slightly farther from the Galactic plane than 
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a: 58 microlensed bulge dwarfs 

b: 166 red giants in BW 
prob=0A9, D=0.16 

c: 220 red clump giants in BW 
prob=0.04, D=0.20 

d: 1813 ARGOS red clump giants 
prob=0. 00, D=0.25 

e: 363 red giants at 6=-li 
prob=0. 00. D=0.33 
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Fig. 13. Cumulative histograms and the Kolmogorov-Smirnow 
tests between the microlensed dwarf sample and the giant star 
samples shown in Fig. [12] The solid thicker red line shows the 
microl ensed bulge dwarfs, (b) the RG B sample from Hill et al.l 
(1201 lh . (c) the red clump sample from lHill etaD (1201 lh . (d) the 
red clump s ample from ARGOS (Ne ss et al.), and (e) the RGB 
sample from lUttenthaler et al.l d2012l) . 



the other samples, this trend could be explained if there were a 
metal-rich population that drops off very rapidly with distance 
fro m the plane. The MDF of 363 RGB stars at (/, b) = (0,-10) 
bv lUttenthaler et ail d2012l) in Fig.fTZe shows a more dominant 
metal-poor peak, possibly confirming the drop-off with \b\ for 
the metal-rich peak. However, contrary to the ARGOS MDF, the 
b = -10° MDF bv lUttenthaler et al.l (1201 2l) contains a larger 
fraction of metal-rich stars ([Fe/H] > +0.3). Since it is located 
even farther from the plane than the ARGOS fields it could mean 
that the ARGOS study suffers from a systematic shift in [Fe/H] 
(i.e., it should be moved to higher [Fe/H]). 

Re gard ing the regions clo ser to the plane, both lRamfrez et all 
(120001) and lRich et al.l (1201 2l) used infrared spectroscopy to study 
M giants in several fields in the innermost 600 pc between 
Baade's window and ~ 150pc (i.e., 4° to 1°) below the plane. 
Neither of these two studies find a major vertical abundance or 
composition gradient. Hence, if there is a gradient, or gradual 
change of stellar populations, it appe ars to disappear cl ose to the 
plane. What is intriguing about the iRich et al] (1201 2l) study is 
that the M giant MDFs peak at [Fe/H] a -0.2 at all latitudes, 
and that there are essentially no M giants with [Fe/H] > 0. This 
is in stark contrast to the findings from both the microlensed 
bulge dwarf sample a s well as t he RG B and RC samples in 
Baade's window from IHill et all (1201 ll) . which all see a large 
fraction of stars at super-solar metallicities. The reason might be 
that very metal-rich M stars actually never reach th e RGB phase 
due to strong stellar winds (see discussion in, e.g.. lCohen et al.l 
I2008I or that the metallicities of M giants may be systemati- 
cally underestimated. In any case, the metal-rich stars are lack- 
ing from studies that probe the MDF using M giants 

To further investigate the variation of the MDF we split 
the microlensed dwarf sample in latitude and then in longitude. 
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(a) 29 microlensed bulge dwarfs 
farther from the plane, |fa|>3' 




(b) 29 microlensed bulge dwarfs 
closer to the plane, |b|^3° 
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(c) 31 microlensed bulge dwarfs 
farther from the centre, M|>2° 




(d) 27 microlensed bulge dwarfs 
closer to the centre, |/|^2° 




Fig. 14. (a) and (b) show the MDFs when the microlensed dwarf 
sample is split into stars farther from the plane (\b\ > 3), and 
closer to the plane (\b\ < 3), respectively, (c) - (d) the MDFs 
when sample is split into stars farther from the centre (|/| > 2), 
and closer to the centre (\l\ < 2), respectively. The curved lines 
represent generalised histograms. 



Figures fT4h and b show the MDFs when the microlensed dwarf 
sample is divided into two samples: one with stars closer than 3° 
to the plane, and one with stars farther than 3° from the plane. A 
two-sided KS-test (prob = 0.32, and D = 0.24) shows that the 
hypothesis that the two subsamples have been drawn from the 
same underlying distribution cannot be rejected. However, to the 
eye, two distributions appear to have distinctly different features, 
and not just a systematic shift which would as in the case of a 
vertical metallicity gradient. The outer sample appears bi-modal 
with a paucity of stars around solar metallicities. The inner sam- 
ple on the other hand contain many stars around solar values. 
Both MDFs span the same metallicity range, but with different 
peaks, indicating a gradual change with b of the fraction of the 
different stellar popul ation components (see also discussion in 
iBabusiaux et al.ll20ld) . When going towards the plane, the frac- 
tion of metal-poor and metal-rich stars drops off while the frac- 
tion of solar metallicity stars increases. The fraction of metal- 
poor stars appears to drop faster than the fraction of metal-rich 
stars. A tentative interpretation could be that there are three pop- 
ulations with different scale heights; could it be the thin disk, the 
thick disk, and a bar population? This interpretation would be in 
line with the results by Ness et al. (in prep.) of their multicompo- 
nent MDF (see Fig.fHb). They interpret their [Fe/H] = +0.14 
peak as belonging to the thin disk surrounding the bulge, the 
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Fig. 15. Age versus [Fe/H] for the microlensed dwarf sample. 



[Fe/H] = -0.23 peak as the true boxy/peanut bulge, and the 
[Fe/H] = -0.60 peak as the old thick disk, maybe being part of 
the bulge. 

Figure [T4"e-f shows the MDFs for the microlensed dwarf 
stars split in two subsamples: one with \l\ > 2 and one with 
\l\ < 2, respectively. The two-sided KS-test yields a similar re- 
sult (prob = 0.32, and D = 0.24) as when splitting the sample in 
latitude, i.e. it is not possible by statistics to reject the claim that 
the two subsamples have been drawn from the same underlying 
distribution. 



4.4. Stellar ages 

Figure[T5]shows the age-metallicity diagram for the microlensed 
bulge dwarfs. The metal-poor stars below [Fe/H] a -0.4 are 
consistently old with ages around 10 to 12 Gyr. The metal -rich 
stars on the other hand show a large dispersion, with ages rang- 
ing from a few billion years up to ~ 13 Gyr. 

The varying age distribution is further highlighted in Fig. IT6l 
which shows the summed age distribution functions (ADF) from 
the individual stars in different metallicity bins (an ADF for 
an individual star include all possible ages from the isochrones 
that are encompassed by the uncertainties of the stellar pa- 
rameters). The summed ADF of the most metal-rich popula- 
tion ([Fe/H] > 0.3) is dominated by young stars and peaks 
around 3-4 Gyr with a long tail towards higher ages. The ADF 
of the most metal-poor stars ([Fe/H] < -0.5) is on the other 
hand dominated by stars older that 10- 12 Gyr. The ADF for the 
stars with -0.5 < [Fe/H] < shows a mixed age structure 
with a broad peak around 5-9 Gyr. The ADF for the stars with 
-0.5 < [Fe/H] < divides into two peaks, one coinciding with 
the metal-rich ADF (~ 4 Gyr) and one coinciding with the metal- 
poor ADF ( 1 0- 1 2 Gyr) . 

4.5. On the presence of young stars in the bulge 

The microlensed dwarf sample signals the existence of a signif- 
icant fraction of low- and intermediate age stars in the bulge. 
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Fig. 17. Y 2 isochrone plots for different metallicities. Four different ages are plotted(l, 5, 10, and 15 Gyr). In each plot, all 58 stars 
in microlensed dwarf sample are included, and those with [Fe/H] within +0.25 dex of the isochrone metallicity have been marked 
by bigger red circles. 
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(a) 14 stars with [Fe/H]>-0.3 
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Fig. 16. Sums of individual age probability distribution functions 
for the microlensed dwarf sample for four metallicity bins (as 
indicated). 



As deep colour magn i tude diagrams of the Galactic bulge (e.g 



Holtzmanetai1ll993 



2003; Clarksonet 



Ortolani et al.l 1 1 995b IFeltzinp & Gilmore 



muT 



Brown et al 



2000; Zoc cali et al 

20101) show a faint red turn-off, indicative of an exclusively old 



and metal- rich population , this is a surprising result. For in- 
stance, the IClarkson et al.l d201 lh study of blue stragglers con- 
cluded that no more than 3 % of the bulge population could be 
younger than 5 Gyr. In our sample of microlensed dwarf stars 
13 out of 58 stars (23 %) have ages lower than 5 Gyr. However, 
if the claimed age uncertainties are taken into account, we only 
have 3 stars out of 58 (~ 5 %) that within 1 <x have an age lower 
than 5 Gyr, which is more o r less in agreement with the estimate 
from IClarkson et al] (1201 lh . It should be noted that our sample 
does indeed contain 13 stars out of 58 (23 %) that have an up- 
per age limit below 9 Gyr, i.e. taking the age uncertainties into 
account these stars are for certain younger than 9 Gyr, pointing 
to a significant intermediate-age population of stars in the bulge. 



While it is feasible to hide a small number of metal-rich and 
young stars "inside" a red and old turn-off, but this cannot be 
large. 

In Fig. [17] we show the microlensed dwarf stars plotted on 
top of isochrones with five different metallicities (as indicated 
on the top of the figure). Each plot contains all 58 microlensed 
dwarf stars, with the stars that have [Fe/H] within +0.25 dex of 
the isochrone metallicity highlighted in red. It is evident that a 
single set of isochrones cannot be used for all stars. Especially 
at low metallicity, many of the stars that have higher metal- 
licities fall outside the isochrones. If we were to assign a sin- 
gle set of isochrones to the sample, say the [Fe/H] = -0.2 
set of isochrones, it would be difficult to claim a significant 
intermediate-age population in the bulge. As we have no stars 
in the upper left corner of the CMD, i.e., within the 1 to 5 Gyr 
isochrones, one would be tempted to claim that the bulge turn- 
off is around 10 Gyr and that the bulge population is all old. This 
is not the case, as can be seen in panel with the [Fe/H] = +0.3, 
many stars fall on young isochrones, stars that would be classi- 
fied as old on a set of more metal-poor isochrones. This illus- 
trates the importance of metallicity as a key ingredient in the 
determination of the age of the bulge. In many cases, such as 
in photometric studies, the metallicities of the stars are lacking, 
and a metallicity has to be assumed. The microlensed stars show 
that even though it is very difficult to claim stars younger than 
say 3-4 Gyr, it is clear that there is a substantial fraction of inter- 
mediate age (5 to 8 Gyr). 

iNataf et a"fl d201 ll) found evidence from studying the red gi- 
ant branch bump that th e bulge should be enr iched in helium. 
This led to an attempt bv lNataf & Gould! (1201 2l) to reconcile the 
age situation by suggesting that the reason for the apparent miss- 
match in age between spectroscopic and photometric results can 
be related to an enhanced He/metals content, Y/Z, in the bulge as 
opposed to that in the Sun. As the stellar isochrones used to in- 
fer the ages are calculated using the standard He abundance the 
resulting spectroscopic ages might be erroneously too low. They 
discuss two systematic effects that arise as a consequence of an 
age offset that should be a function of the evolutionary state of 
the star. The first prediction was that this age offset, when ap- 
plying standard log g — T e g isochrones, should be maximised on 
the subgiant branch and minimised around the main-sequence 
turn-off. The second prediction was that the difference between 
the true absolute magnitude and the spectroscopically inferred 
absolute magnitude should correlate positively with the spectro- 
scopically inferred stellar mass. Using the, at the time, 26 avail- 
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Fig. 18. (a) and (b) show the spectroscopic ages as a function 
of spectroscopic surface gravity for the metal-poor and metal- 
rich bulge dwarfs, respectively, (c) and (d) show the difference 
between absolute magnitude inferred from microlensing tech- 
niques and from spectroscopy as a function of inferred spectro- 
scopic mass for the metal-poor and metal-rich bulge dwarf s, re- 
spectively. (Compare Figs. 2 and 3 in lNataf & Gouldf2012l) . 



able microlensed bulge dwarfs in iBensbv et al.l (1201 Obi 1201 ll) 
they found a positive correlation between the difference in abso- 
lute mag nitude and stell a r mas s for the stars with [Fe/H] > 0. 
This led iNataf & Gouldl d2012l) to tentatively claim that the Y 
abundance should be higher in the bulge than in the Sun. 

The plots in Fig. [T9|are reproductio ns of the bottom plots in 
Figs. 2 and 3 of Nataf & Gouldl (1201 2l) but with the now much 
expanded sample of 58 microlensed bulge dwarfs. Fig.lTSh and 
b show the stellar age versus log g for the metal-poor and the 
metal-rich bulge dwarfs, respectively. Fig. [T8b and d show the 
difference between the absolute / magnitude inferred from mi- 
crolensing techniques and spectroscopy versus the spectroscop- 
ically inferred stellar mass for the metal-poor and metal-rich 
bulge dwarfs, respectively. It is clear that the trend with stel- 
lar mass that was found for the metal-rich bulge dwarfs is no 
longer present. Also, for the stars with sub-solar [Fe/H], where 
they found a negative trend with stellar mass, there is now a pos- 
itive trend. If anything, this shows that a statistically significant 
sample of microlensed bulge dwarf stars is needed to firmly es- 
tablish whether there are any trends or not. 

Figure[l9ji shows the lNataf & Gouldl (1201 2l) prediction (their 
Table 1) of the ratio between the spectroscopically inferred age 
to a "true" age versus surface gravity in the case of a true age 
of lOGyr. As can be seen it would be possible to derive a too 
low He abundance in the models responsible for some of the 
young stars we see in the bulge. For instance, old ages could be 
recovered for some of the stars with ages less the 5 Gyr that have 
log g < 3.8 or \ogg > 4.5 if the He abundance were increased 
by 0.1 or more. However, this is a dramatic change, and there 
would still be stars with ages around, or less than, 5 Gyr between 
4 < logg < 4.3, where an increase in the He abundance will not 
have an equally large impact on the spectroscopically inferred 
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Fig. 19. (a) shows th e predicted ratios (taken from Table 1 in 
INataf & Gouldl d2012l) for a true age of 10 Gyr) between the in- 
ferred spectroscopic age and the true age versus surface gravity 
for four different increases of the He abundance (AY, as indi- 
cated on the right-hand side), (b) shows the age distribution for 
the microlensed dwarf sample, (c)-(e) show how the age distri- 
bution changes with AY if the ages are corrected following the 
relationships in (a). 



stellar ages. In Fig.[T9b we show the age distribution for the 58 
micro lensed dwarf stars and Figs. [T9b-e then show how this 
age distribution changes when the effects predicted in Fig. [T9h 
are applied to the stars. Even though there are stars that become 
substantially older, it is not enough, even in the extreme case 
of AT = +0.098, to to contradict the conclusion that the age 
distribution has a substantial fraction of young and intermediate 
age stars. 

There is other evidence, from AGB star s, of an inter- 
media te age population in the b ulge (e.g., van Loon et al.l 
2003t 



age population in the b ulge (e.g., 
Cole & Weinberg! 120021: lUttenthaler et al l 



2007h . 



Uttenfhaler et al. ( 20071) find evidence for Tc in a sub-sample 
of their C-stars, indicative of third dredge up and a minimum 
stellar mass of 1 .5 M H which implies an upper age limit of 3 Gyr. 
Ivan Loon et al.l (120031) used near infrared CMDs from DENIS 
and ISOGAL to simultaneously derive the extinction, metallicity 
and age of individual stars. They found that the inner region (the 
inner 10°) are dominated by an old population (> 7 Gyr), but that 
an intermediate age population is also present. This is consistent 
with their finding of a few hundred AGB stars with heavy mass 
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Fig. 20. [Fe/H] versus A max for the microlensed dwarf sample. 



loss. Using data from 2MASS, ICole & Weinberg d2002h were 
able to identify a population of carbon stars in the Galaxy that 
traced the bar. Their a nalysis finds that these stars v e ry like ly are 
of intermediate age. iGroenewegen & Blommaertl (120051) also 
found Mira variables with ages of only 1-3 Gyr at latitudes be- 
tween -1.2° to -5.8°. That stars have formed recently "around " 
the bar is not unexpecte d (e.g.. ISanchez-Blazquez et ail 1201 ll) . 
ICole & Weinberg! ((20025) discuss the possibility that the carbon 
stars have wandered in to the bulge/bar region, but are unable to 
draw firm conclusions. 



4.6. The A max puzzle 



Using at the time 16 available microlensing events ICohen et all 
d2010l) reported an unexpected correlation between the metal- 
licity of the microlensed star and the maximum magnification 
of the microlensing event in the sense that metal-rich stars had 
higher maximum magnifications than metal-poor ones. Several 
possibilities were investigated but none could satisfactory re- 
solve the situation. 

Figure l20b shows [Fe/H] versus A max for the now much 
larger sample of microlensed dwa rf stars. Compared to when 
the analysis by ICohen et al.l (1201 Oh was done, the sample now 
contains high-magnification events at low metallicities and low- 
magnification events at high metallicities. However, looking at 
the histograms in Fig. l20k and b, it is clear that the metal- 
rich part of the sample contains a higher fraction of high- 
magnification events than the metal-poor part. Classifying a mi- 
crolensing event with A max > 70 as a high-magnification event, 
and otherwise a low-magnification event (vertical dotted line in 
Fig. l20l . the metal-rich sample contain 84 % high-mag events, 



compared to 51 % for the metal-poor stars. As it is highly un- 
likely that the maximum magnification should have any depen- 
dence on the chemical composition of the background star, this 
correlation must be a result of some selection bias or being due 
to that the environment of the metal-rich stars is different from 
the environment of the metal-poor stars. In the following section 
we will investigate how big such a bias might be and whether the 
differences in the A max distributions between the metal-rich and 
metal-poor samples can be accounted for. 



4.7. Sampling bias? 

The selection of candidate microlensed bulge dwarfs for spec- 
troscopy is made from the OGLE and MOA surveys, which im- 
age through / and R/I filters, respectively. The cutoff in bright- 
ness for effective high dispersion spectroscopy even with 8- 10 m 
telescopes, given that target of opportunity observations cannot 
exceed 1 or 2 hours, is / brighter than about 15.0 mag. 

We attempt to model our sampling process to look for se- 
lection effects. We assume a distance to the bulge of 8.2 kpc 
and a mean extinction of At ~ 1 .0. We adopt v ersion 2 of the 
Yale-Yonsei evolution ary tracks and isochrones (lYi et al.l|2.Q03b 
iDemarque et al. 2004) and select five representative choices of 
[Fe/H], -0.90, -0.43, +0.05, +0.39, and +0.60 dex. Stars with 
subsolar metallicities are assumed to have [a/Fe] = +0.3, while 
higher metallicity stars have no a-enhancement. We adopt the 
Salpeter IMF. 

These isochrones are normalized to have 1000 stars in their 
initial mass function with masses between 0.5 and l.OM . The 
summed stellar mass for each our five metallicity choices for 
isochrones with age 4 Gyr was compared; there is a total range 
of 2 %, which we can ignore for present purposes. 

If we assume that all the bulge stars have the same age, ir- 
respective of metallicity, then at ages older than ~ 5 Gyr, the 
main sequence turnoff is ~0.4 mag fainter for bulge dwarfs with 
[Fe/H] super-solar than it is for dwarfs with -0.9 < [Fe/H] < 
-0.5. This would suggest that in the microlensed sample, the 
intrinsically fainter high-metallicity dwarfs would be under- 
represented. 

However, a key result of this paper is the apparent depen- 
dence of mean age of the bulge dwarfs with metallicity (see 
Fig. [T3t . If we adopt this as true, and allow the stellar mean 
age to decrease as the metallicity increases, with the most metal- 
poor stars having an age of 12 Gyr and the most metal -rich 
ones 4 Gyr, a different result emerges; the younger population 
of super-solar metal-rich bulge dwarfs then has a main sequence 
turnoff that is essentially identical to that of the more metal- 
poor, older bulge populations, with a substantial blue hook as 
is shown in Fig|2T] This is manifested in the star counts in the 
regime 2.2 < M\ < 3.8, which are dominated by the youngest, 
most metal-rich bulge turnoff region stars. In this regime, the star 
counts from this young population are a factor of ~ 2 or more 
higher than those from populations with ages of 10 to 12 Gyr 
(see Fig.l22li. 

When we simulate lensing of this sample, this difference of 
a factor of two is smoothed out somewhat, but is still present 
(see Fig.l23l. At Mi(lensed) = -0.6, corresponding to the limit- 
ing magnitude of our spectroscopic sample, the youngest, most 
metal-rich population has star counts which are about 50% 
higher than those of any of the older, more metal-poor popu- 
lations with the same total stellar mass. Hence we would over- 
estimate the fractional contribution to the total bulge population 
of young, metal-rich stars. For a system of uniform age 10 Gyr, 
the super-solar stars would be under-represented in the lensed 
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Fig. 21. Y 2 isochrones for five different metallicities with differ- 
ent mean ages that decrease with increasing metallicity. 
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Fig. 22. Number of stars along the isochrones shown in Fig. |2T| 



star counts by about a factor of 30 %, and would then be more 
difficult to detect in our microlensed bulge dwarf sample. 

Thus, assuming the adopted metallicity - age relation shown 
in Fig. [15] is valid, the fraction of the bulge stellar population 
which is metal-rich and young is over-estimated in our current 
microlensed bulge dwarf sample of 58 stars by ~50 %. Note that 
a turnoff star from the lower metallicity populations with Mi m 
+3.8 has to be lensed by a factor of 55 to be bright enough to be 
included in the spectroscopic microlensed bulge stellar sample. 

Next, we simulate the distribution of A max in an attempt 
to determine whether selection effects arising from the age- 
metallicity relation of Fig. [15] could be producing the apparent 
deficiency of low metallicity, high A max stars seen in Fig. [20] 

Our simulation for A max is shown in Fig. [24] We find that 
there is a strong trend for the young, metal-rich population to 
dominate over the more metal-poor and older populations which 
is strongest at the lowest magnifications and decreases towards 
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Fig. 24. Simulated magnification distribution for the isochrones 
shown in Fig. [2T1 



higher magnification. Since there are more intrinsically lumi- 
nous stars in the former, they will dominate at the lower magni- 
fication levels. At high magnification, the main sequence turnoff 
can be boosted into the observed sample for the entire range 
of metallicities probed, and hence the bias favoring the high- 
est metallicity young population is reduced. We thus deduce that 
our sample of microlensed bulge dwarfs is biased in favor of 
high-metallicity young stars, and hence our derived MDF is also 
biased. However, it does not appear possible to reproduce the 
apparent relative absence of low-metallicity high magnification 
events with our simulation. 



4.8. Abundance trends 

4.8.1. General appearance for 12 elements 

The plots in Figs [25] and [26] show the abundance trends for O, 
Na, Mg,Al, Si, Ca, Ti, Cr, Ni, Fe, Zn, Y, and Ba for the 58 mi- 
crolensed bulge dwarfs. In Fig.[25]Fe is the reference element 
and in Fig.[26]Mg is the reference element. In all plots the sizes 
of the circles have been scaled with the ages of the stars (larger 
circles equal higher ages). The a-element plots in Fig. [25] (i.e., 
[O, Mg, Si, Ca, Ti] versus [Fe/H]) show a plateau of elevated 
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Fig. 27. Abundance trends for the cr-elements Mg, Si, Ca, and Ti. Left-hand-side panels show [X/Fe] versus [Fe/H], and right-hand 
panels show [Fe/X] versus [X/H]. Filled grey circ les mark the micro l ensed bulge dwarfs and the Red and blue circles are nearby 
thick and thin disk dwarf stars, respectively, from dBensbv et al.ll2003l I2005L and in prep.). 



[c/Fe] abundance ratios at metallicities below [Fe/H] « -0.4, 
signalling fast chemical enrichment by massive stars. These stars 
are all old. The a-trends then decline towards solar values, and 
even though most of the stars are old, there are one or two that 
start to show up with lower ages. At super-solar metallicities the 
trends generally level out, meaning that there is some kind of 
equilibrium between SN II and SN la enrichment. The age struc- 
ture is also very complex with a mix of stars at all ages, 2 to 
15Gyr. What is notable is that the abundance trends are very 
well-defined and generally have very low scatter. 

By using an a-element that mainly comes from a single 
source (in this case core-collapse supernovae), one might get 
a clearer picture on the chemical history of a stellar popula- 



tion. The "x-axis" could then be interpreted as a proxy for time. 
Hence, Fig.|26]shows the abundance trends again, but now with 
Mg as the reference element. Especially interesting plots are now 
the Fe, Cr, and Ni plots where we can see flat, under-abundant, 
[(Fe, Ni, Cr)/Mg] trends for low [Mg/H] values. This is con- 
sistent with rapid enrichment by massive stars. The [(Fe, Ni, 
Cr)/Mg] ratios start to increase around [Mg/H] « -0.2 to —0.1, 
signalling that low-mass stars are becoming the dominant source 
for chemical enrichment of the interstellar medium. 
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4.8.2. Fast enrichment - revising the position of the knee 



In iBensbv et al.l (1201 Obi) we noted that even though the metal- 
poor bulge dwarf stars generally follow the abundance trends 
outlined by the thick disk it was also apparent that they were 
slightly more a-enhanced, mainly being located on the upper 
rim of the thick disk abundance trends. Figure [27] shows the 
abundance trends for the tr-elemen ts Mg, Si, Ca , and Ti with 
the nearby thick disk stars from (IBensbv et alJ I2003L I2005L 
an d in prep.) included for comparison purposes. As found 
by IBensbv & Feltzin ill (120121) . discriminating between thin and 
thick disk stars based on stellar age appears to produce "cleaner" 
results when applied to local samples, and we have therefore 
marked stars with ages greater than 9 Gyr (and age uncertainties 
less than 4 Gyr) as thick disk stars (red circles), and stars with 
ages less than 7 Gyr (and age uncertainties less than 4 Gyr) as 
thin disk stars (blue dots). 

At low [Fe/H], the metal-poor bulge dwarfs show an a- to- 
iron ratio that is similar to what is seen the nearby thick disk 
(whi ch also is similar to what is observed in the stellar halo, see, 
e.g.. lNissen & Schustertl20Tol) . However, it is possible that the 
declining [or/Fe] ratios (or rising [Fe/or] ratios) occur at slightly 
higher [Fe/H] (or [a/H]) than for the local thick disk. As we now 
have many more microlensed dwarf stars around and slightly 
below solar [Fe/H], where the thick disk [a/Fe] trends decline, 
we are more confident to claim that this slight shift might be real. 
This slight shift in the position of the "knee" is clearly evident 
in Mg and Ti, slightly less in Si, and barely in Ca. If it is real it 
can be interpreted as being due to the bulge having experienced 
a faster enrichment than the local thick disk, leading to the onset 
of enrichment from low-mass stars occur at higher metallicities. 

Does this mean that the possible connection between the 
bulge and the thick disk becomes weaker? As of now we are 
only able to compare the bulge with the thick disk in the solar 
neighbourhood while it would be desirable to compare with the 
thick disk in the inn er regions of the Galax y. The only such study 
to date is the one bv lBensbv et al.l d2010al) but their sample of 44 
red giants is unfortunately too small to whether or not the inner 
thick disk is different from the nearby thick disk. If there is a ra- 
dial metallicity gradient in the thick disk it would mean that the 
inner thick disk regions have experienced a more rapid evolu- 
tion, and hence pushed the "knee" to higher metallicities. What 
we see in the bulge is a "knee" that starts to decline perhaps 
0. 1 dex higher in [Fe/H] than seen in the nearby thick disk. 



5. Discussion 

The first microlensed dwarf stars that were studied provided a 
real surprise, an MDF heavily skewed to super-solar metallici- 
ties an d significantly different from what had b een found from 
giants dJohnson et al.ll2007l: ICohen et al.ll2008l) . The difference 
was somewhat elevated by the very first of our VLT ToO obser- 
vations from which we found many of the microlensed dwarfs to 
have super-solar metallicities. But as the number of microlens- 
ing events grew we started to also get metal-poor d warfs in larger 
numbers, and a clearly bi-modal MDF emerged (IBensbv et alJ 
1201 Obi 1201 ll) . That MDF was still very different from the best 
MDF from high- resolution spectra o f a representative sample 
of bulge giants dZoccali et al.l l2008|). A recent re -analysis of 
the IZoccali et al.l d2008l) data by iHill et ail (1201 ll) and further 
increased number statistics of the microlensed dwarf stars has 
changed this, and there is now good agreement between the over- 
all appearance of the MDFs traced by microlensed dwarfs and 
the giant stars in the Galactic bulge (see Figs. IT21 and IT31). 
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Fig. 28. [Ba/Mg] versus [Fe/H] for the micr olensed dwarf sam- 
ple. T he two-component bulge model by iTsuiimoto & Bekkl 
(|2012|) is shown (blue dashed line is the first component that 
forms quickly, and red solid line is the second component that 
forms from pre-enriched gas). 



The MDF is an importa nt constraint to models of galac- 
tic chemical evolution (e.g.. iMatteucci & Francois] 1 19891) . The 
MDF of the Galactic bulge is best probed through the detailed 
elemental abundances of stars. Most studies use the intrinsi- 
cally bright red giant stars for this. The first studies showed a 
metal-rich population with a somewhat broadened distribution 
of [Fe/H]. Coupled with the h igh a abundances this pointed 
to a short formation time-scale (Rich 119881: iMcWilliam & Rich] 
ll994HFulbright et al.ll2007h . However, the mean [Fe/H] changed 
between these studies leaving substantial room for interpreta- 
tion. Rece ntly, with the arrival of the multi-fibre spectrograph 
FLAMES dPasquini et al.l2002l) on the VLT, it has become feasi- 
ble to probe further down the red giant branch and reach warmer 
stars that are less pro ne to evolutionary effects (compare lack of 
metal-rich M giants in lRich et al.l2012|) and also to the red clump 
which should give a reliable tracer of the MDF. 

Chemical evolution models have been created to attempt to 
explain the broad spread in [Fe/H] and to explain the fact that 
the bulge appears to have two populations: one with a low mean 
[Fe/H] and one with super-solar [Fe/H]. Fo r example, the two- 
compo nent chemical evolution model by ITsuiimoto & Bekkl 
(|2012|) is able to explain the [Ba/Mg] versus [Fe/H] pattern seen 
in the dwarf stars (Fig. 1281 as well as a bi-modal MDF. They 
also find that a top-heavy IMF is essential to explain the high 
me tallicity of the metal -rich component of the MDF. 

iGrieco et al.l d2012l) also explores a multi-component model 
of the chemical evolution of the Galactic bulge. They model a 
two component bulge with one representing the classical bulge 
and the other the pseudo-bulge (bar population?). One prediction 
from their model is that there could be an abundance gradient 
within the classical, spherical bulge but that also the differences 
between the chemical evolution of the spherical component and 
the bar component could be (erroneously) interpreted as a gradi- 
ent. 

In simulations of the formation and evolution of the Milky 
Way in a cosmological context, the star formation of the 



bulge lasts for longer than a sing l e starburst ( Rahimi et all 
| 20 1 Ob iKobavashi & Nakasatol 1201 ll IScannapieco 



etal. 2011 



iDomenech -Moral et al. 2012 ). For example, w hile most of the 
bulge stars in IKobavashi & Nakasatol d2011l) are old (80%> 
10 Gyr, and 90% > 8 Gyr), 10% of stars form later, at super-solar 
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metallicities. While [Kobavashi & Nakasatol (201 1|) simu l ated a 
sing le history for a MW-like galax y, IScannapieco et all (1201 ll) 
and Domenech-Mor al et alj (l20T2h each modelled several halo 
histories. The galaxies showed different star formation histories 
for all components, and it is therefore likely that star formation 
histories are not universal in each galaxy type. However, sev- 
eral of their simulationsshow extended star formation histories 
for the bulge, although the interpretation of IScannapieco etall 
d201 lb is complicated by the fact that their "inner spheroid" ex- 
tends ~ lOkpc. While it is not yet clear how well these simu- 
lations model the behaviour of baryons through their sub-grid 
physics, that some fraction of the bulge consists of stars younger 
than 8 Gyr is not surprising within the framework of current the- 
ories. 



So what is the bulge? Where do all the observational evidence 
point to? One thing is clear, there is overwhelming evidence for 
more than one population in the bulge. Whether it is two, or even 
more, it is a fact that the inner region of the Galaxy is where the 
major Galactic stellar populations (if they all exist in this re- 
gion) should meet and overlap. For instance, the interpretation 
by Ness et al. (in prep.) of their multicomponent MDF is that the 
[Fe/H] = +0.14 peak belongs to the thin disk surrounding the 
bulge, the [Fe/H] = -0.23 peak is the true boxy/peanut bulge, 
and the [Fe/H] = -0.60 peak is the old thick disk, which may or 
may not be part of the bulge. Adding the Galactic bar (possibly 
formed from disk stars), and the influence that it has on its sur- 
roundings, the situation is truly intricate. However, at the same 
time the metal-poor bulge shows abundance trends that are very 
well-defined with low scatter, which could argue against a bulge 
composed of more than two populations (as proposed by Ness et 
al. in prep., see Fig.[T2tl). 

As noted in Sect. 14.8.21 the or-element trends for the metal- 
poor bulge dwarfs now appear to be slightly different from 
those of the nearby thick disk, possibly implying that the metal- 
poor bulge population experienced a somewhat more intense 
star formation than the local thick disk. This does not necessar- 
ily weaken a possible connection between the metal-poor bulge 
component and the thick disk. It might just instead be a reflec- 
tion of the more intense star formation in the inner regions of the 
Galaxy than in and around the solar circle. 

The abundance trends seen for the metal-rich bulge is 
partly consistent with what is seen for metal-rich nearby thin 
disk dwarf stars (e.g.. feensbv et all I2003L 12005b iReddv etalJ 
l2006t IBensbv et al1l2007bh . although disk stars with metallici- 
ties higher than [Fe/H] > +0.4 are rarely seen. Given the wide 
range in ages, the metal-rich bulge could be the most complex 
region in the Galaxy, possibly hosting stars from several popula- 
tions. 

A bimodal age distribution in the bulges of external spi- 
rals galaxies with bars has been recently found through a stellar 
population synthesis analysis bv lCoelho & Gadottil (1201 ll) . with 
mean ages of 4.7 and 10.4 Gyr. Interestingly, unbarred galaxies 
do not show this bimodal distribution, perhaps implying that bars 
have the effect of rejuvenating bulges. The above results are in 
line with our findings based on Galactic bulge dwarf stars, and a 
larger sample will allow us to obtain firmer results for compari- 
son with external bulges. 



6. Summary 

In this study we have presented a detailed elemental abundance 
analysis of 32 microlensed dwarf and subgiant stars in the bulge. 



Together wi th the previous sample of 26 microlensed bulge 
dwarfs from IBensbv et al.l (1201 Obi 1201 ll) the sample now con- 
tains 58 stars. In summary, the findings and main results from 
this data set are: 

- The metallicities span the full range between [Fe/H] = -2 
to +0.6. The MDF that appeared clearly bi-modal, even with 
a paucity of stars around [Fe/ H] = 0, when based on the 
26 stars in IBensbv et al.l d20 lib , is now less so. The overall 
shape of the MD F is now s i milar to the MDF in Baade's 
window found by iHill et all d201 lb . as traced by both red 
giants and red clump giants. Furthermore, signatures of more 
than two components is starting to emerge in the microlensed 
dwarf star MDF. 

- The microlensed dwarf stars more metal-poor than [Fe/H] < 
-0.4 are all old with ages around or greater than 10 Gyr. The 
more metal-rich stars, on the other hand, have a wide age 
distribution with a significant fraction of stars having ages 
lower than 5 to 6 Gyr, and quite a large number of stars with 
intermediate ages around 5-7 Gyr. These are essentially in- 
distinguishable from an old turn-off when the metallicities 
are unknown. We have investigated in some detail whethe r 
an He enriched bulge, as proposed bv lNataf & Gouldld2012l) . 
could explain the large fraction of young stars in the mi- 
crolensed dwarf sample, but the effects are not large enough. 
We have also investigated sampling biases and find that our 
microlensed dwarf sample could be biased in the sense that 
the number of young metal-rich stars is over-estimated (per- 
haps by as much as 50 %). In any case, neither sampling bias, 
nor a He enriched bulge, can fully account for the large frac- 
tion of young and intermediate age stars in the microlensed 
bulge dwarf star sample. 

- With the now much expanded sample of 58 stars we see 
that the metal-poor bulge population have abundance trends 
slightly different from what is observed in the nearby 
Galactic thick disk. The "knee" in the [a/Fe] - [Fe/H] trends 
occurs at a slightly higher metallicity in the bulge, indicat- 
ing that the metal-poor bulge component has experienced a 
somewhat faster star formation rate. 

- The microlensed bulge dwarfs with solar and super-solar 
[Fe/H] show abundance trends similar to what is seen in the 
nearby metal-rich thin and thick disks. If the bar population 
is formed from disk material, this is what could be expected. 

Even though it is difficult to pinpoint the exact origin of the 
bulge, observational evidence now clearly indicates that the 
bulge has a complex structure with wide age and metallicity dis- 
tributions. Combined with the cylindrical rotation of the bulge, 
the bulge appears to be a conglomerate of Galactic stellar popu- 
lations under the influence of the bar. 
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